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ABSTRACT 


Pert'nrmance  data  from  the  NASA-Army  OH-6A  higher  harmonic  control 
(HHC)  flight  test  program  showed  significant  reductions  in  main  rotor  shaft  torque 
and  engine  power  in  hover  and  forward  fiight  flj  The  unsteadv  aerodynamics  with 
the  higher  harmonic  control  system  application,  including  wake  elTecis.  were  con¬ 


sidered  to  study  whether  such  power  reductions  are  feasible 

An  airfoil  oscillating  in  pure  plunge  can  achieve  propulsive  Ibrcc  ("  Kalzmayr 
effe:-."),  as  in  the  case  of  birds’  wings  flapping  in  flight.  Here  it  will  be  shown 
tha  this  effect  can  be  enhanced  in  the  presence  of  layers  of  shed  vorticiiy  with 
th,'  proper  phasing,  In  addition,  while  it  is  known  that  an  airfoil  oscillating  in 
pitch,  will  typically  produce  drag  at  most  values  of  reduced  frequency,  it  is  found 
that  the  presence  of  another  layer  of  shed  vorticity  of  the  proper  phase,  can 
reduce  the  drag  on  the  pitching  airfoil  depending  upon  wake  spacing,  reduced  fre- 


qUciiCy,  and  j;iia5c. 


V 


iiuci  buiiic  cuiiuiiunib  me  auueu  layci  ui 


I avers  ot 


SllCli  VlJl- 


ticity  will  even  result  in  propulsive  force  acting  on  the  pitching  airfoil  similar  to 
the  "Kalzmayr"  effect  for  the  plunging  case. 

It  was  found,  for  the  OH-6A  helicopter,  that  the  measured  reductions  in 
main  rotor  shaft  torque  and  engine  power  are  feasible  when  evaluated  with  respect 
to  the  "Katzmayr  effect"  and  the  additional  drag  reduction  or  propulsive  force 
obtained  due  to  pitch  and  plunge  oscillations  with  the  effect  of  adjacent  wake 
layers  of  shed  vorticity. 
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Summary 

open  loop  performance  data  from  the  NASA-Army  OH-6A  higher  harmonic 
control  (HHC)  flight  test  program.  Wood  et  al.  [I],  show  significant  reductions 
in  main  rotor  shaft  torque  and  engine  power  in  the  airspeed  regime  from  hover 
to  100  knots.  Depending  upon  HHC  controller  phase  and  helicopter  airspeed, 
reductions  in  power  were  recorded  as  large  as  20%,  with  reductions  of  the  order 
of  10%  being  typical. 

Higher  harmonic  control  is  an  active  helicopter  vibration  control  concept 
which  alters  the  aerodynamic  loads  on  the  rotor  loades  such  that  the  blade  re¬ 
sponse  is  reduced.  This  in  turn  reduces  the  vibratory  forces  and  moments  acting 
at  the  hub,  which  cause  airframe  vibration.  Basically,  HIiC  is  an  electronic, 
computer-controlled  active  vibration  suppression  system  which  senses  and  cancels 
vibrations  in  a  helicopter  by  N/revolution  feathering  or  pitch  motion  of  the  rotor 
blades,  N  being  the  number  of  blades’.  Verv  substantia!  vibration  reduction  was 
achieved  under  a  NASA-Army  sponsored  program  using  a  motdfied  OH-6A  he¬ 
licopter.  Higher  harmonic  blade  pitch  control  was  implemented  by  superimposing 
4/rev.  (32  Hz)  swashplate  motion  upon  basic  cyclic  and  C(  .lev  live  control  inputs, 


'  li.Kobcns  Wood.  R,  W  I’owcrs.  J  1 1.  Clinc.  and  C  U.  1  Irrninoiid.  i  >ii  I'-'cw’'. oping  and 
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During  the  program,  in  addition  to  reducing  vibration  levels,  the  HHC  system 
showed  the  potential  for  decreased  helicopter  power  requirements.  Recent  study 
of  the  OH-6A  flight  test  data  indicates  that  the  mechanism  by  which  the  power 
reduction  was  achieved  is  related  to  the  unsteady  aerodynamics  associated  with 
HHC.  This  is  supported  by  the  fact  that  higher  harmonic  control  by  the  very 
nature  of  the  method  by  which  it  is  implemented  through  oscillating  an  airtbil 
in  pitch  with  resulting  plunge  motion,  requires  definition  of  unsteady  aerodynam¬ 
ics  to  adequately  model  its  physics’. 

The  fundamental  closed  form  solutions  of  Theodorsen  [3],  and  Loewy  [4] 
provide  the  basis  for  theoretical  work  in  this  area.  The  closed  form  theory 
shows  rapid  changes  in  the  lift  deficiency  function  with  changes  in  reduced  fre¬ 
quency,  wake  spacing  and  frequency  ratio.  In  the  past,  emphasis  in  the  study 
of  unsteady  aerodynamics  has  tended  to  focus  on  flutter  instability  and  the  effea 
of  unsteady  aerodynamics  on  generating  lift  and  torsional  loads. 

In  this  research,  we  are  interested  in  performance,  and  emphasis  will  be  on 
the  effect  of  unsteady  aerodynamics  on  the  drag  of  the  airfoil.  When  the  drag 
force  is  reversed,  acting  to  propel  the  airfoil  forward,  it  is  generally  referred  to 
as  either  negative  drag  or  propulsive  force.  The  classic  reference  on  this  subject 
is  that  by  Garrick  [5] 


'  h.  R.  Wood,  Max  !■  Plat/cr.  Ahmed  Abourahma.  Mark  A.  Couch  "  Cn  llic  Unsteady 
Aerodynamics  of  Iligtier  llamiome  Control."  Paper  No  C  17.  Nineteenth  Puropean 
Rotoreralt  I'orum  Ceniobbio  tConio).  Ita!\,  September  1 99.1. 
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While  Garrick’s  [5]  work  shows  that  an  airfoil  oscillatinu  in  pitch  will 
typically  nroduce  drau  in  the  lower  reduced  frequency  range,  k.  it  is  found  that 
the  presence  of  another  layer  of  shed  vorticity  of  the  proper  phase  can  reduce 
the  drag  on  pitching  airfoils  and  depending  upon  wake  spacing,  reduced  frequen¬ 
cy,  and  phase,  may  even  enhance  the  propulsive  force  acting  on  the  airfoil  similar 
to  the  "Katzmayr "  effect  for  plunging  airfoils. 

The  specific  mechanisms,  that  have  been  used  in  this  research  to  discuss 
'he  phenomenon  of  OH-6/\  helicopter  rotor  power  reduction,  are:  ( 1 )  the  unsteady 
-odynamics  associated  with  a  plunging  airfoil  and  how  that  produces  a  propul¬ 
sive  force  or  decrease  in  rotor  power;  (2)  the  intluerce  of  a  layer  of  wake 
vorticity  on  this  propulsive  force,  which  represents  the  influence  of  the  unsteady 
aerodynamics  of  one  rotor  blade  on  adjacent  rotor  blades,  (3)  the  unsteady  aero¬ 
dynamics  associated  with  an  airfoil  oscillating  in  pure  pitch  and  how  this  pro¬ 
duces  in  general  a  drag  force  on  the  airfoil  unless  accompanied  by  layers  of  shed 
wake  vorticity  ;  and  (4)  the  influence  of  the  phasing  of  a  layer  of  shed  vorticity 
on  this  drag  force,  and  how,  depending  upon  specific  phasing,  the  effect  on  the 
airfoii  will  be  either  drag  or  propulsive  force. 

The  analysis  of  the  above  effects  was  based  upon;  (1)  the  classical  flat 
plate  unsteady  aerodynamic  theory  for  an  incompressible  fluid  of  Theodorsen  [3], 
its  extension  to  include  chordwise  forces  (drag)  by  Garrick  [5],  the  rotor  wake 
analysis  of  Loewy  [4],  (2)  the  panel  method  developed  by  Plat.c.'  et  al.  [6], 
that  can  be  applied  fo’"  unsteady  incompressible  How  past  airfoils  or  airfoil  com¬ 
binations  of  arbitrary  geometry  and  amplitudes  of  different  kinds  of  motions. 
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In  the  first  approach.  Theodorsen  s  well  known  lift  deticiency  lunction  C(k) 
has  been  modified  by  Loevw  tor  the  case  of  a  hovering  rotor.  Ai’doil  drag  or 
propulsive  force  was  determined  by  the  methods  of  Garrick  [2],  In  the  second 
approach,  the  panel  code  allowed  a  systematic  study  of  the  lift  and  drag  produced 
bv  oscillating  airfoils  and  airfoil  combinations  in  an  incompressible  invisctd 
flow. 

The  panel  code  results  were  checked  against  different  references.  The  results 
showed  '.rood  agreement  with  these  references.  For  the  parameters  ot  the  OH- 
6.A.  rotor  in  hover,  the  main  conclusion  was  that  the  reductions  in  drag  .and 
propulsive  forces  generated  on  an  airfoil,  oscillating  in  pitch  or  plunge,  due  to 
wake  phasing,  are  sufficient  to  explain  the  reductions  in  power  recorded  on  the 
OH-6A  during  the  N.ASA-Army  flight  test  program  [1] 


I .  INTRODUCTION 


A.  GENERAL 

The  determination  of  the  unsteady  aerodynamics  associated  with  a  helicopter  rotor  has 
presented  a  challenge  to  the  helicopter  analyst  since  the  conception  of  the  helicopter.  The 
effect  of  blade  oscillations  on  resulting  airloads  has  prompted  new  interest  in  oscillatory 
blade  aerodynamics  in  connection  with  possible  performance  gains.  Rotar\'  wing  unsteady 
aerodynamics  is  considerably  more  complex  than  fixed  wing  aerodynamics. 

The  simplest  case  is  that  of  hover,  but  even  here  the  problem  is  complicated  by  the 
influence  of  the  shed  helicoptei  wake  from  the  rotor.  In  forward  flight,  the  problem  is 
funher  complicated  by  the  fact  that  the  velocity  of  a  local  blade  element  ranges  at  the 
blade  tip  from  transonic  on  the  advancing  side  to  low  subsonic  on  the  retreating  side,  and 
at  the  blade  root  from  low  subsonic  on  the  advancing  side  to  reversed  flow  on  the  retreat¬ 
ing  side. 

For  this  reason  there  is  a  great  need  for  advances  in  the  present  theory  in  order  to 
explain  modem  helicopter  problems,  yet,  the  number  of  advances  in  rotorcraft  unsteady 
aerodynamic  theory  are  very  limited,  so  we  need  to  review  current  unsteady  aerodynamic 
theory  for  rotorcraft  and  to  advance  the  theory  in  one  of  the  areas  where  improvement  is 
badly  needed. 

The  goal  of  this  thesis  is  to  explain  the  measured  power  reductions  in  the  OH-6A 
higher  harmonic  control  (IDIC)  flight  test  program  conducted  by  C  R.  Wood  et  al  [1]. 
This  research  will  consider  the  unsteady  aerodynamics  with  HTIC  system  application, 
including  wake  effects,  in  order  to  determine  whether  such  power  reductions  are  feasible. 
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B.  ROTARY  WING  UNSTEADY  AERODYNAMICS 

In  this  section,  some  elements  of  unsteady  airfoil  aerodynamics  will  be  discussed. 
Consider  an  airfoil  undergoing  sinusoidal  pitch  and  plunge  motion.  As  the  airfoil  pert'orms 
an  oscillation,  vortices  are  shed  into  the  medium  with  a  circulation  strength  equal  in  mag¬ 
nitude  to  the  increase  in  circulation  about  the  airfoil,  but  opposite  in  direction.  These  dis¬ 
turbances  are  stored  in  the  lluid  because  the  shed  vorticity  convccts  downstream  at  the 
local  flow  field  velocity 

The  counter-rotating  ^ortices  induce  a  sinusoidal  llow  field  which  further  changes  the 
.1.  The  result  of  this  complex  flow  field  is  a  time  difference  or  delay  between  the 
airfoil’s  motion  and  the  induced  aerodynamic  forces.  This  delay  is  known  as  the  phase 
lag,  VI/  [34] .  To  simplify  calculations  for  this  type  of  motion,  it  is  common  to  describe 
the  airfoil  position  and  the  associated  aerodynamic  forces  by  complex  variables.  For  pure 
plunge  oscillations  the  vertical  motion  of  the  airfoil  is  described  by  the  real  part  of  the 
following  equation; 

hit)  =  h„  e'"‘  . (1.1) 

where  1\,  is  a  complex  number,  and  o  is  the  frequency  of  oscillation.  Similarly,  the  lift  is 
Hc’"'  d  by  ; 


L=  L^,  r  e 


1  (Dt 


where  U,  is  the  quasi-steady  lift  given  by  the  expression; 


(1.2) 


L,  =  .5  pU“S  (  h  /  U) 


(1.3) 
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This  is  termed  quasi-steady  because  the  angle  of  attack  is  represented  by  h/U. 
The  values  of  q,  and  \|/  represent  the  magnitude  and  the  phase  respectively,  of  the 
true  instantaneous  lift  reladve  to  the  quasi-steady  lift.  The  variables  r,  and  \i/  in 
general  depend  on  the  reduced  frequency  k,  the  Mach  number  M.  and  the  Reynolds 
number  For  inviscid,  incompressible  fluid  the  values  of  i;,  and  \\i  will  be  only  a 
function  of  k  [34],  .A  complete  solution  for  the  oscillating  ftat-plate  airfoil  in  an 
incompressible  inviscid  flow  has  been  obtained  by  Kussner  and  Theodorsen  and 
reproduced  from  Reference  [35]  in  Figure  (1.1  ). 
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C.  ADDITION  OF  WAKE  EFFECT 

For  the  helicopter  in  hover,  the  rotation  of  the  wing  introduces  a  number  of  features 
that  require  special  attention,  notably  the  returning  vortex  wake,  the  time  vaiying  free 
stream,  and  radial  flow,  and  a  fundamentally  transcendental  geometry  that  requires  either 
approximate  or  numerical  solutions.  The  lift  on  a  wing  is  due  to  its  bound  circulation. 
Conservation  of  vorticity  in  the  flow  requires  that  there  be  trailed  and  shed  vorticity  in 
the  wake  behind  the  wing. 

The  spanwise  variation  of  bound  circulation  results  in  trailed  vorticity  parallel  to  the 
free  stream  direction.  Time  variation  of  the  bound  circulation  leads  to  shed  vorticity  par¬ 
allel  to  the  wing  span.  The  wake  is  composed  of  sheets  of  vorticity  convected  downstream 
from  the  trailing  edge  by  the  free  stream  velocity  The  wake  of  a  rotor  in  hover  or  in  ver¬ 
tical  flight  consists  of  helical  vortex  sheets  below  the  disk,  one  fro^i  each  blade.  Unsteady 
motion  of  the  rotor  blade  will  produce  a  shed  vorticity  in  the  wake  spirals.  With  low  disk 
loading  { at  low  collective  pitch  settings ),  the  wake  remains  near  the  rotor  disk  and  there¬ 
fore  passes  close  to  the  following  blades  Thus  the  wake  vorticity  is  not  convected  down¬ 
stream  of  the  airfoil  as  with  a  fixed  wing,  and  the  shed  vorticity  sheets  below  the  rotor 
disk  must  be  accounted  for  to  correctly  estimate  the  unsteady  loads. 

For  high  inflow  (high  collective  pitch  settings  )  or  in  high-speed  forward  flight,  the 
wake  will  be  convected  away  from  the  blades.  The  returning  shed  wake  influence  is  pri¬ 
marily  a  concern  of  hover  and  vertical  flight,  and  forward  flight  at  airspeeds  less  than 
100  knots.  Assuming  high  aspect  ratio  of  the  blade,  lifting  line  theory  requires  a  knowl¬ 
edge  of  the  loads  on  the  blade  section,  and  the  returning  shed  wake  of  the  rotor  must  be 
incorporated  into  the  two-dimensional  unsteady  airfoil  theory. 
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The  wake  far  from  ihe  blade  section  will  have  little  influence,  so  the  emphasis  is  on 
modeling  the  wake  near  the  blade,  which  for  low  inflow  consists  of  vortex  sheets  that  are 
nearly  planar  surfaces  parallel  to  the  disk  plane.  Based  un  these  considerations,  a  two- 
dimensional  model  for  the  unsteady  aerodynamics  of  the  rotor  can  be  constructed.  Loewv 
in  1957,  developed  a  two-dimensional  model  for  the  unsteady  aerodynamics  of  the  blade 
in  a  hovering  rotor.,  including  the  effect  of  the  helical  shed  wake. 

D.  HIGHER  HARMONIC  CONTROL  (HHC) 

The  primary  source  of  helicopter  vibration  is  the  higher  harmonic  blade  loads  gener¬ 
ated  by  the  aerodynamic  environment  at  the  rotor  disk.  For  an  N-bladed  rotor,  the  oscil¬ 
latory  (N-l)P,  NP,  and  (N+1)P  blade  loads  are  transmitted  as  exciting  hub  forces  and 
moments  to  the  airframe  at  a  dominant  frequency  of  N  per  rev  (NP).  A  significant  reduc¬ 
tion  in  vibration  level  has  been  achieved  by  industry  in  the  last  twenty  five  years.  The 
desired  goal  was  to  achieve  0.02  g  and  this  goal  cannot  be  met  w’  hout  a  quantun^ 
advance  in  vibration  control  technology.  The  vibration  levels  first  specified  by  the  U.  S 
Army  for  the  AAHAJTTAS  procurement  in  the  early  1970s  were  later  revised  upward  in 
the  mid  1970s  to  better  reflect  realistic  design  goals  consistent  with  the  state-of-the-art  in 
helicopter  vibration  control  [1]. 

A  higher  harmonic  control  system  drives  the  blades  in  pitch  at  the  (N-l)P,  NP,  and 
(N+1)P  harmonics  of  rotor  rotational  speed,  generating  new  unsteady  airloads,  which  in 
combination  with  the  new  oscillatory  inertial  loads,  cancel  the  harmonics  of  blade  loads 
that  cause  airframe  vibration.  Therefore  the  vibrations  are  suppressed  at  the  source. 
Although  there  are  several  ways  to  implement  HHC  on  a  rotor  system,  the  approach 
generally  followed  to  date  has  been  by  blade  root  feathering  using  swashplate  oscillations . 
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By  means  of  electro-hydraulic  serv'o-actuarors.  the  swashplate  is  excited  in  the 
collective,  longitudinal  cyclic,  and  lateral  cyclic  modes  at  NP  resulting  in  blade  pitch 
oscillations  at  three  distinct  frequencies  of  (N-l)P.  ,MP,  and  (N+1)P  in  the  rotating  frame. 
In  the  system  considered,  developed  for  the  four-bladed  OH-6A  helicopter,  higher  har¬ 
monic  blade  feathering  for  vibration  reduction  is  achieved  by  superimposing  4/rev  swash- 
plate  motion  upon  basic  collective  and  cyclic  Bight  control  inputs. 

Perturbing  the  stationary  swashplate  at  4/rev  both  vertically  and  in  pitch  and  roll 
results  in  third,  founh.  and  fifth  harmotiic  blade  feathering  in  the  rotating  system.  Fourth 
harmonic  blade  feathering  is  achieved  by  oscillating  the  swashplate  vertically  about  its 
collective  position,  while  third  and  fifth  harmonic  blade  feathering  in  the  rotating  system 
results  from  4/rev  tilting  of  the  stationary  swashplate  in  pitch  and  roll  about  its  cyclic  tilt 
position  [2].  Results  of  experimental  efforts  by  Hammond  [52]  showed  that  successful 
suppression  of  vibration  can  be  achieved  by  osciilating  the  blades  at  relatively  small 
angles 

E.  OH-6A  (NASA/ARMY/HUGHES)  TEST  PROGRAM  DATA 

The  goal  of  this  research  has  been  a  systematic  study  of  the  OH  6A  helicopter 
HHC  flight  test  program,  in  order  to  understand  and  explain  the  significant  reduction  in 
main  rotor  shaft  torque  and  engine  power  that  was  measured  in  fight.  Shown  in  Figure 
(1.2)  are  the  power  reduction  results  from  the  HHC  Open  Loop  Flight  Test  Program 
(NAS/^J  Army/  Hughes  -  Sept.  1982). 
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figure  (  1-2  ):  Powei*  KcUiirtioii  Kesuits  IVoiii  lliglier  llaiiiKiiiic  C'oiilrol  (IIIIC) 


Open  Loop  I’esi  I’rogiani  (NASA/Ariiiy/lliiglics  -  Sepi,  r>82) 


l  igiiics  (13)  and  (14)  (rioiii  Kef.  [ I ])  show  the  very  substantial  vibration  reduction 
achieved  under  a  NASA-Anny  sponsored  program  using  a  modified  OH-6A  helicopter.  Higher 
hartnoiiic  blade  pitch  control  was  achieved  by  superimposing  4/rcv,  (32  Hz)  swashplale 
motion  upon  basic  cyclic  and  collective  control  inputs.  The  aircraft  was  flown  from  zero 
airspeed  to  1(10  knots  with  the  HHC  system  operated  both  open  loop  (manually)  and  closed 
loop  (computer  controlled). 


Ainr.rf  Fti  ik  iasi 


I'igiirc  (  1.3  ):  (iomitarison  of  4/icv  Vertical  Vibration  at  (lie  Pilot 
Seat  Wrsiis  Airspeed  with  IIIIC  "on"  and  "ofT'. 
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Figure  ( 


1.4):  Comparison  of  4/rcv.  Lateral  Vibration  at  Die  Pilot  Seat 
Versus  Airspeed  with  IIHC  "on"  and  "ofr'. 
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The  significance  of  power  reduction  means  that:  (I)  the  helicopter  payload  car  be 
increased.  (2)  the  maneuverability  can  be  markedly  increased.  (3)  marginal  operating  coiidi- 
tions,  such  as  at  sea  in  high  sea  states  can  be  better  met.  and  (4)  the  helicopter  operating 
costs  can  be  substantially  reduced  Increased  interest  in  HHC  and  its  potential  benefits  Tor 
the  helicopter  now  require  that  we  try  to  use  the  proper  tool  to  understand  this  phenome¬ 
non.  It  would  appear  that  we  presently  have  the  computational  tools,  namely  the  paiv'l 
method  for  unsteady  inviscid,  incompressible  ilow  past  airfoils  or  airfoil  combinations  of 
arbitrary  geometry,  developed  by  Platzer  ct  al  [6],  and  the  single  wake  closed  fonn  anal¬ 
ysis  developed  from  Loewy’s  theory  by  Wood  and  Couch  [40], 
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F.  SCOPE 


Chapter  II  contains  a  quick  review  of  the  most  important  research  starting  with  the 
work  of  Knoller[10]  in  1909  and  Betz  [II]  in  1912.  They  explained  the  bird’s  ability  to 
generate  forward  thrust  by  means  of  wing  flapping.  In  1922.  R.  Katzmayr,  experimentally 
verified  the  work  of  Knoller  and  Betz.  Birnbaum,  von  Karman.  Theodorsen.  and  Garrick  in 
the  years  1924  through  1936,  were  the  pioneers  of  flapping  wing  aerodynamics. 

The  fundamental  closed  form  solutions  by  Theodorsen  [3],  Loevvy  [4],  Garrick  [5],  pro¬ 
vide  our  basis  for  theoretical  work  in  this  field.  While  Loewy's  work  on  wake-induced 
flutter  helps  to  explain  the  phenomenon,  it  also  points  up  difficulties  to  be  overcome.  That 
is  the  closed  form  theory  shows  rapid  changes  in  the  lift  deficiency  function  with  changes 
in  reduced  frequency,  k;  wake  spacing,  h;  and  frequency  ratio,  m. 

Garrick  further  developed  the  work  of  Theodorsen  by  deriving  the  horizontal  force  for¬ 
mulas  for  airfoils  or  airfoil-aileron  combinations  oscillating  in  any  of  the  three  degrees  of 
freedom,  vertical  flapping,  torsional  oscillation  about  a  fixed  axis,  and  angular  oscillation 
of  the  aileron  about  a  hinge.  Also  presented  in  Chapter  II,  Loewy's  two-dimensional  model 
for  the  unsteady  aerodynamics  of  the  blade  of  a  hovering  rotor,  is  presented  to  proviae  the 
base  for  discussing  the  obtained  results.  Also  presented  in  Chapter  II.  is  a  quick  survey  of 
the  most  relevant  references  who  developed  and  used  panel  methods  to  study  unsteady 
aerodynamics  of  oscillating  airfoils  in  an  incompre,.uble.  inviscid  flow. 
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In  the  third  chapter,  the  propulsive  force  is  estimated  according  to  the  OH''6A  helicop¬ 
ter  parameters  obtained  from  reference  [36],  In  this  analysis,  the  blade  is  divided  into  sec¬ 
tions.  Assuming  hovering  conditions,  the  reduced  frequency  at  each  section  was  evaluated. 
Using  the  unsteady  panel  code  the  time  history  of  the  airfoil  drag  (  or  propulsive  force  ) 
coefficient  was  recorded.  For  the  assumed  plunging  motion  with  plunging  amplitude  as 
obtained  from  reference  [  1  ],  for  a  given  mode  of  vibration,  the  average  propulsive  force 
for  each  section  was  calculated. 

The  total  average  propulsive  foice  was  obtained  by  adding  together  the  contribution 
from  each  blade  section.  The  obtained  results  showed  that  the  resulting  propulsive  force 
from  a  single  airfoil  oscillating  in  plunge  is  not  enough  to  explain  the  significant  power 
reduction  on  the  OH-6A  helicopter.  This  led  us  to  smdy  the  aerodynamic  interference 
between  two  blades,  which  is  the  main  subject  of  the  fourth  chapter. 

In  Chapter  IV,  tile  two-airfoil  unsteady  panel  code  developed  by  Pang  [  25  ],  is  used 
to  study  the  wake  interference  between  two  airfoils  This  code  has  been  formulated  to 
solve  for  the  potential  flow  for  two  airfoils  executing  unsteady  motion  in  an  inviscid 
incompressible  flow  medium..  It  is  an  e.xtension  of  the  code,  developed  by  Teng  [24],  i'or 
single  airfoils  The  technique  uses  the  well  known  panel  method  for  steady  How  and 
extends  it  to  unsteady  flow  by  accounting  for  the  continuous  shedding  of  vortices  into  the 
trailing  wake. 
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Numerous  case-runs  are  presented  to  illi:sirate  the  aerodynamic  interterence  between  two 
oscillating  airtbils.  Results  obtained  for  p'unging  and  pitching  airfoils  show  a  significant 
production  of  propulsive  force  at  certain  wake  phase  angles.  In  Chapter  V’  the  question  of 
propulsive  efficiency  due  to  plunging  or  pitching  and  the  possibility  of  single-degree-of- 
freedom  pitching  flutter  are  discussed.  Chapter  VI  discusses  the  findings  of  this  research 
and  presents  the  main  conclusions. 
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11.  BACKGROUND 


The  eftect  of  plunge  and  pitch  oscillations  on  the  resultant  drag  of  an  airfoil 
has  been  the  object  of  aeroelastic  research  back  almost  to  the  time  of  the  Wright 
brothers  (1903).  In  1909  and  1912.  respectively,  Knoller  flO],  and  Betz  [11]  were 
the  first  ones  to  explain  the  birds’  ability  to  generate  a  forward  thrust  by  means 
of  wing  flapping’ .  .\nother  one  of  the  earliest  investigators  of  this  era  was  R. 
Katzmayr  [7],  who  in  1922  conducted  two  senes  of  experiments  to  explore  the 
eiTects  of  oscillatory  fluid  or  airfoil  motion  on  resultant  lift  and  drag. 

Katzmayr’s  experiments  which  were  performed  in  the  aerodynamic  laboratory 
of  the  Technical  University  of  Vienna,  may  be  divided  into  two  series.  In  the  first 
series,  the  angle  of  attack  of  the  wing  model  v/as  changed  by  oscillating  the  wing 
model  about  an  axis  parallel  to  the  span  and  alright  angles  to  the  airflow.  In  the 
second  series  of  experiments,  the  airfoils  were  stationary  while  the  airflow  itself 
was  subjected  to  periodic  oscillations.  Here,  for  certain  air  velocities  and  frequen¬ 


cies,  and  considering  thick  airfoil  sections.  Katzmayr  measured  negative  drag  or 


positive  propulsive  force.  This  discovery,  to  be  later  known  as  "Katzmayr  effect". 


helped  to  explain  the  propulsive  mechanism  of  birds  flying  through  the  air. 


'  I’lal/.cr.  M  i'..  Ncacc.  K  S  .  ;ind  I’aiig,  C  K  ." Aerodynamic  Analysis  ol'  l•lapplne  Wing 
I'ropuision ".  .AIAA  puixn  No  V.t-04X4.  .tlsl  .Aerospace  .Sciences  inccliiig,  l<cno.  NV,  .lanuaiy 
11-14.  199.1 
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Later.  E.  G.  Richardson  [8]  used  the  "Katzmayr  erfect"  to  explain  the  locomotion 
of  fish  through  water.  In  his  paper,  Katzmayr  summarized  his  results  to  show:(l) 
the  effect  of  flowing  air,  whose  direction  is  undergoing  constant  periodical  change, 
is  quite  favorable,  (2)  the  wing  section  which  exhibits  favorable  characteristics  in  a 
constant  airflow,  works  better  in  an  oscillating  flow. 

In  1935,  Theodore  Theodorsen  [3]  introduced  the  general  theoiy  of  aerody¬ 
namic  instability  and  the  mechanism  of  flutter.  In  the  first  part  of  his  paper. 
Theodorsen  developed  the  velocity  potentials  due  to  flow  around  oscillating  air- 
foil-ailerons.  In  the  later  part  of  the  paper,  he  developed  the  differential  equations 
of  motions  and  determined  the  flutter  speed. 

The  most  relevant  work  with  respect  to  the  present  research  is  that  by  I.  E. 
Gamck  [5],  who  applied  the  work  of  Theodorsen  [3]  to  obtain  a  closed  form  ana¬ 
lytical  expression  for  the  propulsive  force  generated  by  an  oscillating  airfoil.  An 
airfoil  oscillating  in  pure  plunge  or  flapping  motion  was  shown  by  Garrick  to  have 
a  propulsive  force  throughout  the  entire  range  of  reduced,  k,  while  for  an  airfoil 
oscillating  in  pitch,  propulsive  force  can  only  be  achieved  for  a  specified  range  of 
reduced  frequency,  k,  where  k  must  be  greater  than  2.0  . 
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In  1942,  VV,  Schmidt  [1 4]  recognized  that  the  Knoller-Betz  or  Katzmayr  effect 
applies  to  both  a  flapping  airfoil  in  uniform  flow  or  a  stationary  airfoil  in  oscillat¬ 
ing  flow.  He  demonstrated  experimentally  that  a  stationary  airfoil  positioned  in  the 
wake  of  a  flapping  airfoil  increases  the  propulsive  efficiency  to  almost  100  percent 
throughout  the  whole  frequency  range  because  the  stationary  airfoil  converts  the 
vortical  energy  generated  by  the  flapping  airfoil  into  additional  thrust.  He  also 
found  that  the  mechanical  limitations  inherent  in  pure  flapping  motions  could  be 
overcome  by  an  arrangement  which  he  called  the  wave  propeller. 

H.  Bosch  [15]  invest' A  the  interference  of  two  lifting  surfaces  in  two- 
dimensional  unsteady  ucompressible  flow,  using  classical  flat  -plate  theory.  He 
showed  that  a  harmonically  flapping  airfoil  upstream  of  a  stationary  airfoil  increases 
the  propulsive  efficiency  to  almost  100  %,  thus  confirming  Schmidt’s  experimental 
findings.  Bosch’s  analysis  is  limited  to  flat-plate  airfoils  oscillating  about  a  zero 
incidence  mean  position. 

In  1957,  Loewy  [4]  developed  a  two-dimensional  model  for  the  unsteady 
aerodynamics  of  the  blade  of  a  hovering  rotor.  He  derived  an  expression  for  the 
two-dimensional  oscillatory  lift  deficiency  function  as  a  function  of  the  reduced 
frequency,  wake  spacing  anu  the  ratio  of  the  airfoil  frequency  of  oscillation  and  the 
rotational  speed. 
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The  wake  in  hover  or  vertical  flight  consists  of  helical  vortex  sheets  below  the 
disk,  one  from  each  blade.  Unsteady  motion  of  the  rotor  blade  will  produce  shed 
vorticity  in  the  wake  spirals. 

In  Figures  (2.1  A)  through  Figure  (2. ID),  from  Ref  [4],  Loewy  illustrated  the 
basic  elements  of  the  unsteady  flow  field  of  a  helicopter  rotor.  rotor  in  hover¬ 
ing  or  vertical  (light  produces  a  trailing  tip  vortex  with  a  downward  a.xial  veloc¬ 
ity  that,  if  otherwise  undisturbed,  forms  a  contracting  helix  as  shown  in  Figure 
(2.1  A).  If  the  inflow  over  the  rotor  disk  is  constant,  then  the  fluid  off  the  trailing 
edge  of  the  blades  makes  a  helical  surface  with  horizontal  radial  elements  as  shown 
in  Figure  (2. IB). 

If  there  is  an  oscillation  in  blade  effective  angle  of  attack,  blade  lift  will  alter¬ 
nate  also,  and  as  a  result  of  these  changes  in  lift,  vortices  will  be  shed  continu¬ 
ously  at  the  blade  trailing  edge.  These  vortices  fall  along  the  horizontal  radial 
elements  of  the  helical  surface  shown  in  Figure  (2. IB),  so  long  as  the  oscillations 
in  angle  of  attack  are  small.  Figure  (2.iC)  illustrates  this  local  sheet  of  shed  vor¬ 
ticity,  It  should  be  noted  that  vorticity  is  considered  to  be  on  the  helical  surface 
shown  in  Figure  (2. IB);  the  vertical  displacements  from  that  surface  shown  in  Fig¬ 
ure  (2. 1C)  represent  the  strength  of  the  vorticity  at  a  particular  azimutltal  and  radial 
position  [4]. 
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The  variation  in  this  venical  displacement  (hence  vortex  strength)  around  the 
azimuth  corresponds  to  the  history  of  the  motion  of  a  given  blade  element  at  a 
fixed  radius.  Variation  of  shed  vonex  strength  in  the  radial  direction  at  any  fixed 
azimuth  angle  is  a  function  of  the  variation  with  blade  span  of  (1)  blade  planform; 
(2)  amplitude  of  oscillation  of  effective  angle  of  attack,  and  (3)  relative  air  veloc¬ 
ity.  Since  (from  the  Helmholtz  theorem  )  vorticity  cannot  begin  or  end  in  space,  a 
variation  of  shed  vorticity  in  the  radial  direction  implies  that  there  are  trailing  vor¬ 
tices  at  constant  radii  similar  to  and  inboard  of  the  lip  vortex.  These  trailing  vorti¬ 
ces  have  been  induced  in  Figure <2. 1 D). 
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The  schematic  drawings,  Figures  (2.1  A)  through  Figure (2. ID),  Ref.  [4],  indicate 
the  difficulty  of  attempting  to  obtain  a  complete  representation  of  unsteady  rotor 
aerodynamics.  One  means  of  simplifying  the  problem  is  to  consider  extreme  "pitch" 
values  for  the  helix.  For  the  case  of  ve.ry  high  inflow  velocity,  u,  in  relation  to  Q, 
and  the  opposite,  where  u  is  very  low  compared  to  £1  This  is  shown  in  Figure 
(2.2). 

When  the  vertical  spacing  between  adjacent  helicai  surfaces  of  shed  vorticity  is 
very  large,  then  we  expect  that  all  shed  vorticity  beyond  a  small  fraction  of  a  rev¬ 
olution  would  be  too  far  below  the  blade  in  question  to  have  a  significant  effect. 
Under  these  conditions,  it  would  be  sufficient  to  account  for  only  the  attached 
vortex  sheet  withi.n  that  fraction  of  a  revolution,  as  in  Figure  (2.2A), 

On  the  other  hand,  when  the  vertical  spacing  between  the  adjacent  helical 
surfaces  of  shed  vorticity  is  very  small,  all  the  sheets  of  shed  vorticity  tend  to  pile 
up  on  each  other,  and  the  effect  of  that  vorticity  close  to  the  blade  in  question 
(  shed  bv  blade  nasses  )  is  of  more  imnortance  than  that  which  exists  bevond 
a  reference  azimuth  angle  on  either  side  of  the  blade,  This  situation  is  depicted  in 
Figure  (2,2B),  Ref  [4], 
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A.  High  Inilow 


II.  Low  Inilow 


I  igurc  (2.2)  Sclicmalie  llf|ncscii<ji<ion  of  Hiislcady  Rolor  Flow  FipJd. 
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With  low  rotor  thmst  values  {  low  collective  pitch  )  the  wake  remains  near 
the  rotor  disk  and  therefore  passes  cl  .se  to  the  following  blades.  Thus  the  shed 
vorticiiy  sheets  below  the  rotor  disk  tiiu...  be  accounted  for  to  correctly  estimate  the 
unsteady  loads.  For  high  inflow  or  high-speed  forward  flight  the  rotor  wake  will 
be  convected  away  from  the  blades.  The  wake  far  from  the  blade  section  will  have 
little  influence,  so  emphasis  is  on  modeling  the  wake  near  ihe  blade  In  the  case 
of  low  inflow,  this  consists  of  vortex  sheets  that  are  nearly  planar  surt'aces  paral¬ 
lel  to  the  disk  plane  [4] 

Based  on  these  considerations.  two-dimensional  model  for  the  unsteady  aero¬ 
dynamics  of  the  rotor  can  be  constructed.  It  is  assumed  that  the  chord,  amplitude 
of  the  oscillation  in  etfcctive  angle  of  attack,  and  relative  airspeed  vary  slowly 
enough  with  span  so  that  what  occurs  acrodynamically  at  one  blade  radius  station 
is  essentially  duplicated  on  either  side  of  it. 

The  assumption  meanr  that  the  flow  problem  at  a  given  blade  section  is  two- 
dimensional.  Furthermore,  consistent  with  the  idea  that  only  the  vorticity  near  the 
blade  section  has  an  important  clfcct,  one  may  allow  the  planar  rows  of  vorticity  to 
extend  to  infinity  in  the  horizontal  direction  in  order  to  achieve  mathematical  sim¬ 
plification.  It  should  be  stated  that  this  is  a  first-order  theory,  so  that  such  effects 
as  that  of  the  wake  upon  itself  have  been  neglected,  fhe  oscillatory  flow  induced 
at  the  airfoil  in  the  free-steam  direction,  for  example,  will  affect  the  airfoil  lift  and 
moment  as  well  as  the  vorticity  shed  at  the  trailing  edge. 
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Since  the  shed  vorticity  is  being  accounted  for  as  perturbations  in  the  main 
How.  this  component  of  induced  velocity  will  be  small  compared  to  the  free-steam 
velocity  (except  when  the  spacing  approaches  zero)  and  must  react  with  small,  first- 
order  airfoil  displacements.  Because  of  the  physical  unrealities  associated  with  the 
standard  mathematical  model  representation  of  the  problem  (such  as  zero  thickness 
of  wake  and  airfoil,  zero  displacement  of  the  wake  from  horizontal  planes,  and 
infinite  velocities  at  the  core  of  a  vortex),  the  accuracy  of  results  as  the  wake 
spacing  goes  to  zero  is  questionable.  The  effect  of  the  horizontal  component  of 
induced  flow,  therefore,  has  been  neglected  in  Loewy's  analysis. 

Shown  in  Figure  (2.3)  is  the  two-dimensionalized  model  of  unsteady  rotor  flow 
[4],  The  portion  of  the  circular  surface  which  is  determined  by  (1)  a  particular 
blade  radius,  (2)  the  azimuth  angle  on  either  side  of  the  blade  section  (within 
which  the  shed  vorticity  is  of  importance),  and  (3)  the  vertical  distance  spanned  by 
a  given  numbers  of  rows  of  vorticity,  can  be  regarded  as  a  plane,  one  in  which 
the  two-dimensional  unsteady  aerodynamic  problem  may  be  attacked. 

This  is  shown  in  Figure  (2.3  A),  with  an  arrow  indicating  that  the  rows  of 
voiticity  under  the  rotor  disc  extend  to  infinity.  Consistent  with  the  idea  that  only 
the  vorticity  near  the  blade  section  has  an  important  effect,  one  may  allow  the  pla¬ 
nar  rows  of  vorticity  to  extend  to  infinity  in  the  horizontal  direction  in  order  to 
achieve  mathematical  simplification.  These  final  modifications  are  shown  in  Figure 
(2.3  B),  which  is  a  two-dimensional  model.  Ref  [4],  of  unsteady  rotor  aerodynamics 
for  a  single-bladed  rotor  operating  at  low  inflow. 
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The  general  approach  used  by  Loevvy  is  the  same  as  that  of  Schwarz  [30] 
Only  the  effects  resulting  from  oscillatory  motion  of  a  thin  airfoil  are  considered, 
and  the  airfoil  itself  is  represented  by  continuous  vorlicity  distributed  along  a 
straight  line  [9],  It  may  be  also  noted  that  the  effect  of  viscosity  in  damping  the 
shed  vorticity  could  be  accounted  for  by  considering  only  a  finite  number  of  rows 
beneath  the  reference  airtbil,  or  by  multiplying  the  vortex  strength  by  an  experi¬ 
mental  decay  function  related  to  the  time  through  the  rotational  speed  and  the 
number  of  rotational  cycles  corresponding  to  a  particular  vortex  row. 
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This  has  been  omitted  from  the  analysis  since  vortex  decay  time  is  iong 
compared  to  the  rotational  periods  of  most  rotors.  An  extension  to  Loewy ’s  work  is 
that  by  Shipman  and  Wood  [29]  in  1970.  The  authors  present  a  method  for 
determining  rotor  blade  flutter  in  forw'ard  flight.  Incompressible  unsteady  aerody¬ 
namic  theory  was  applied  where  shed  vorticity  both  in  and  below  the  plane  of 
the  rotor  is  accounted  for.  This  was  made  possible  by  noting  that  the  damping  for 
rotor  blade  aeroelastic  instability  varies  with  velocity.  Since  the  tangential  velocity 
of  any  blade  segment  varies  with  azimuth,  so  vdll  flutter  damping'. 

It  is  assumed  at  the  onset  of  flutter  that  oscillations  will  begin  to  build  up 
prior  to  the  blade  reaching  a  critical  azimuth  position,  then  decay  as  the  blade 
moves  beyond  this  point.  This  build  up  and  decay  means  that  the  vorticity  shed 
due  to  the  oscillations  will  be  contained  within  a  double  azimuth  region  on  either 
side  of  the  critical  azimuth  position.  Assuming  this  region  to  be  small  allows  the 
wake  system  to  be  two-dimensional. 

The  lift  deficiency  function  resulting  from  the  two-dimensional  wake  model  is 
compared  with  earlier  results  obtained  for  a  helicopter  in  hover  and  fixed  wings.  In 


the  limiting  case  when  the  advance  ratio  is  very  small,  Shipman  and  Wood’s  liri 
deficiency  reduces  to  Loewy’s  lift  deficiency  function.  Also  in  the  limit  as  the 
wake  spacing  is  very  large,  this  function  reduces  to  Theodorsen’s  lift  deficiency 
function  . 


'  Shipman,  K.  W.  ,  and  I:.  R.  W<kkI  "  A  Two-Dimensional  Thuoiy  for  Rotor  Rladc  Iduttcr 
m  touvurd  I'liplu.  "  .lounial  ol'  Aircral't,  Vol.  No.  S.  Number  12.  pp.  1008-1015. 

December  1971, 
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The  theory  was  applied  to  bending  torsion  lluiter  tor  the  tip  segment  of  a  rotor 
blade.  The  results  followed  the  normal  trends  of  having  the  tlutter  velocity  decrease 
as  the  blade  center  of  gravity  was  moved  aft  and  of  having  the  flutter  velocity 
increase  as  the  blade  stiffness  increased  . 

The  previously  shed  wakes  are  destabilizing  so  that  they  reduce  the  flutter 
velocity.  The  buildup  and  decay  produces  wakes  that  are  essentially  centers  of  vor- 
ticity.  Thus  their  position  with  respect  the  reference  blade  will  be  sensitive  to  the 
advance  ratio  which  determines  the  vertical  spacing.  It  was  found  that  as  the 
wakes  v/ere  brought  closer  by  reducing  the  advance  ratio  and/or  inflow  ratio,  the 
flutter  speed  is  decreased  [29] 

The  phenomenon  of  airfoil  flutter  is  caused  by  energy  extraction  out  of  the 
wind  stream  by  means  of  airfoils  which  are  free  to  oscillate  in  pitch  and  plunge. 
As  described  by  Duncan  [19],  this  effect  can  be  easily  demonstrated  in  a  wind  tun¬ 
nel  if  a  device  is  built  which  allows  to  vary  the  phase  relationship  between  pitch 
and  plunge.  Phase  angles  ne?'"  90  degrees  produce  maximum  energy  extraction  so 
that  the  device  operates  as  aflutter  engine.  McKinney  and  DeLaurier  [20]  demon¬ 
strated  that  this  effect  can  be  used  to  convert  wind  energy  into  mechanical  energy. 
Scherer,  McKinney  and  Delaurier  used  Theodorsen's  flat-plate  theory  for  perfor¬ 
mance  estimates  of  the  flapping  foil  propulsor  and  of  the  wingmill,  respectively. 

It  was  shown  by  W  S»-hmidt,  that  optimum  energy  extraction  occurs  for  har¬ 
monically  Tupping  airfoils  upstream  of  a  stationary  airfoil  lltere  is  a  need  for 
aerodynam’c  analysis  methods  which  enable  systematic  studies  of  the  effect  of  air¬ 
foil  geometry,  incidence  angle  variation,  amplitude  of  oscillation  and  aiifoil  interac¬ 
tion  on  the  aerodynamic  forces  due  to  incompressible  flow  over  airfoils  or  airfoil 
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combinations  executing  pitch  or  plunge  motions.  A  panel  method  was  developed  in 
order  to  do  parametric  studies  of  these  effects.  In  the  past,  a  number  of  investiga¬ 
tors  have  solved  the  steady  How  problem  using  source  and  vortex  paneling,  the 
most  prominent  ones  being  Hess  and  Smith  [21],  .4  few  authors  have  extended  this 
approach  to  the  case  of  unsteady  motion  of  single  airfoils,  notably  Basu  and  Han¬ 
cock  [22]  and  Kim  and  Mook  [23]. 

Basu  and  Hancock  [22]  presented  a  numerical  method  for  the  calculation  of 
the  pressure  distribution,  force,  and  moments  on  a  two-dimensional  airfoil  undergo¬ 
ing  an  arbitrary  unsteady  motion  in  an  inviscid  incompressible  flow.  The  method 
was  applied  to  (i)  a  sudden  change  in  airfoil  incidence,  (ii)  an  airfoil  oscillating  in 
high  frequency  and  (iii)  an  airfoil  passing  through  a  sharp-edged  gust 

Teng  [24],  developed  a  computer  code  for  the  numerical  solution  of  unsteady 
inviscid  incompressible  flov/  over  an  airfoil.  Pang  [25]  extended  the  work  of  Teng 
and  developed  a  computer  code  to  study  unsteady  airfoil  interference  effects.  The 
technique  uses  the  well  known  panel  methods  for  steady  Bow  and  extends  it  to 
unsteady  flow  by  accounting  for  the  continuous  shedding  of  vortices  into  the  trailing 
wake. 
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111.  DRAG  REDUCTION  DUE  TO 


BLADE  FLAPPING 
(KATZMAYR  EFFECT) 

A.  INTRODUCTION 

In  this  chapter,  a  discussion  will  first  be  presented  of  the  physics  of  flapping 
airfoils  which  explains  the  bird’s  ability  to  generate  a  forward  thrust  by  means  of 
wing  flapping.  This  will  be  followed  by  a  description  of  two  methods  capable  of 
predicting  the  pressure  distributions,  lift,  pitching  moment,  and  drag  or  thrust  on  airfoils 
which  execute  harmonic  oscillations  in  plunge  or  pitch. 

The  first  method  is  based  on  the  generalization  of  the  well-known  steady-state 
panel  method  to  the  case  of  unsteady  airfoil  motion.  It  is  based  on  the  assumption 
of  inviscid,  incompressible  flow,  The  second  method  is  based  on  the  classical  The- 
odorsen  theory  for  oscillating  flat-plate  airfoils  in  incompressible  inviscid  flow.  The 
unsteady  panel  method  permits  the  investigation  of  finite-thickness  airfoils  of  arbitrary 
geometry  which  may  execute  rather  general  unsteady  motions  ( ramp  and  oscillatory  ). 
In  contrast,  Theodorsen’s  analysis  is  restricted  to  thin  airfoils  ( flat-plates )  and  to  small 
amplitude  oscillations.  The  chapter  ends  with  a  comparison  of  results  computed  by  both 
methods. 
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B.  THE  PHYSICS  OF  FLAPPING  AIRFOILS  IN  LOW-SPEED  FLIGHT 


The  phenomenon  of  thnjst  generation  due  to  wing  flapping  or  due  to  steady  flight 
in  a  sinusoidal  gust  car.  be  explained  in  a  rather  elementaiy  way.  Consider  the  latter 
case  of  flight  through  a  sinusoidal  gust,  as  shown  in  Figure  (3.1).  The  relative  wind 
vector  changes  direction  between  a  maximum  positive  and  negative  angle  of  attack.  The 
resultant  aerodynamic  lift  is  perpendicular  to  the  instantaneous  wind  vector.  If  we  as¬ 
sume  the  viscous  drag  to  be  small,  then  it  is  easily  seen  that  the  sinusoidal  gust  gen¬ 
erates  a  sinusoidally  varying  thrust. 

A  similar  effect  is  generated  if  the  airtbil  executes  a  sinusoidally  varying  plunge 
motion  (  wing  flapping )  about  an  otherwise  steady  flight  conditions.  This  explanation 
of  thrust  generation  due  to  wing  flapping  or  due  to  flight  through  a  sinusoidal  gust  was 
first  advanced  by  Knollei  [10]  and  Betz  [11].  Its  expeiimental  demonstration  was 
accomplished  by  Katzmayr  [  7  ]  and  is  therefore  generally  referred  to  as  Katzmayr 
effect  ’ . 

Since  it  is  well  known  that  every  change  in  aiifoil  incidence  causes  the  shedding 
of  a  starting  vortex,  a  more  precise  explanation  of  the  Katzmayr  effect  has  to  take  into 
account  the  continuous  vortex  shedding  due  to  wing  flapping  or  due  to  flight  through 
a  sinusoidal  gust. 


’  IMatzcr.  M  !•  ,  Ncucc.  K  S  ,  and  l’an(z,  C  K./'Acrodyrianin;  Analysis  ol'  l  lapinny  Wiiif’  Propulsion", 
jiajic  1 .  AIAA  paper  No  91-04X4,  .list  Aciosjiaec  .Sciences  ineeling,  Reno.  NV.  .lanuaiv  11-14.  199.1 
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Iliriisl 


Figure  (3.1)  Pro|Milsivc  force  on  plunging  airfoil 
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Similarly,  the  quantitative  prediction  of  the  induced  oscillatory  forces  ( lift  and 
thrust )  had  to  await  the  development  of  a  complete  unsteady  aerodynamic  theory. 
Such  a  theory  was  first  accomplished  in  the  1930’ s  by  Theodorsen  [3]  in  the  United 
States  and  by  Kiissner  [16]  in  Germany.  It  can  be  found  in  the  textbooks  on  aeroelas- 
ticity,  for  example  the  book  by  Bisplinghoff,  Ashley,  Halfman  [41].  Therefore  only 
the  final  formulas  will  be  given  here  and  the  reader  is  referred  to  those  texts  for 
additional  details.  More  recently,  a  second  approach  to  the  determination  of  the 
unsteady  aerodynamic  forces  on  oscillating  airfoils  became  feasible  with  the  advent  of 
high  speed  computers  and  the  development  of  efficient  numerical  solution  methods.  The 
fundamental  building  blocks  of  this  method  are  presented  in  the  next  two  sections. 

C.  PANEL  CODE  FOR  STEADY  INCOMPRESSIBLE  INVISCID  FLOW 

In  potential  flow  theory,  the  flow  field  around  an  airfoil  may  be  represented  by 
the  velocity  potential.  Considering  contributions  from  the  free  stream  flow  and  the 
source  and  vorticity  distribution,  the  total  potential  may  be  constructed; 

(p  =  (pp  -  (Ps  r  Ov  (3.1) 

where 


<Pp  =  U  [  X  cos  (  a  )  +  y  sin  ( a  )  ] 


(3.2) 


33 


The  source  distributions  q(s)  vary  i’rom  panel  to  panel,  while  the  vorticity  strength 
y  (s)  is  assumea  constant  for  all  panels.  The  value  of  represenbng  the  How  past  an 
airfoil  by  surt'ace  singularity  distributions  lies  in  the  fact  that  these  singularity  distri¬ 
butions  automatically  satisfy  Laplace's  equation,  the  governing  flow  equation  for  invis- 
cid  incompressible  flow: 


O  -=  d>  =0  (3,3) 

Since  Laplace's  equation  is  a  linear  homogeneous  second  order  parial  differential 
equation,  the  superposition  principle  used  in  Equation  (3.1)  holds.  The  boundary  con- 
dilioiis  include  ilov/  langency  at  control  points  (niid-poinis  of  panels)  and  the  Kutta 
condition  at  the  trailing  edge,  requiring  equal  tangential  velocities  for  the  first  and  last 
panels.  By  evaluating  the  integrals  along  the  airfoil  surface,  the  potential  may  be  de¬ 
termined  at  anv  point  in  the  flow  field.  Each  point  is  defined  at  a  radius  (r)  and  angle 
(0)  from  a  chosen  reference  point  on  the  airfoil. 


T'hp  rpi-£).r<ar</-»fa  tKic  i  c*  flip  liior^in^T  Tho  Qir+V\il  »c 

1  i  iw  I  wiviiW  fuOiiil  111  eiiiij  -j  cliv*  v  ii  i  wu\i*  iiiw  mi  k  i  t  ^ 

by  a  number  of  defined  points,  called  nodes  [34]  More  points  produce  greater  resolution 
and  better  accuracy.  One  hundred  to  two  hundred  points  are  usually  sufficient,  with 
the  larger  number  used  for  more  complicated  airfoil  shapes  or  more  involved  calcula¬ 
tions,  TTie  lines  connecting  these  nodes  are  the  panels.  There  are  n+l  nodes  with  first 
and  last  node  overlapping. 
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I'igiirc  (3.2)  depicts  (he  panel  gcomc!r\'.  Numbering  starts  at  the  trailing  edge,  then 
progresses  along  the  lower  surface.  leading  edge,  and  upper  surface,  and  ends  at  the 
trailing  edge.  The  unit  normal  vectors  (nj)  are  perpendicular  to  the  panels  and  directed 
outward  from  the  airfoil  surface.  The  unit  tangential  vectors  (tj)  are  parallel  to  (he 
panels  and  (he  positive  direction  is  defined  with  increasing  numbering  (n  to  n^-l  ) 

I  he  panels  may  vary  in  length,  with  the  e.xception  of  the  first  and  last  panels,  which 
must  be  equal  in  (^rder  to  use  (he  Kutta  condition  at  the  trailing  edge. 


Figure  3.2  Panel  Method  Geometry 
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The  required  input  consists  of  the  number  of  nodes  on  the  airibii  surtace.  the 
coordinates  of  the  nodes  referenced  from  the  leading  edge,  and  the  angie  of  attack  In 
degrees.  The  program  produces  normalized  \elocities  and  pressure  coefficients  at  each 
control  point  as  output.  The  use  of  intluence  coefficients  ieads  to  a  straightforward 
procedure  for  programming  the  equations. 


An  aerodynamic  intluence  coetficiem  is  defined  as  the  '.eiocitv  inducea  at  a  field 
point  by  a  unit  strength  singularity  distribution  on  one  panel.  For  the  two-dimensional 
steaay  tlow  problem,  t.ne  following  intluence  coefficients  are  needed  1.'41 


.A‘\,  :  normal  velocitv  comoonent  induced  at  the  1  ’’  nanei  control  noint  bv  unit 


strength  source  distribution  on  the  panel 


li  +  l 


A..  =  sin  (0  - 0.) In  --  -  -rcos(0  -0)3- 
1.1  '  1  r  r  '  '  ' 

u 

=  i/2 


1 


:  =  1 


(3,4) 


A^j  ;  tangential  velocity  component  induced  at  the  i'''  panel  control  point  by  unit 
strength  source  distribution  on  the  i'*’  panel 


a'.  =  --  sin  (0.  -  0.)  p..  -  cos  (0. -0  )ln~ 
"  2:t,  ’  h  '  .1  r 


=  0 


'’■1  +  I 


,  I 


(3.5) 


1  =  1 


B",j  .  normal  velocity  component  induced  at  the  i''’  panel  control  point  bv  unit 
strength  vorticity  distribution  on  the  j'*’  panel 

r 

F'TJ  _  sin  (9  -  0.)  t 

u 


B‘‘.  =  -  cos  (0.  -  0.)  In sin  (0  -  0.)  B.. 

11  '»n:*  ‘  ’  r  '  ' 

ij 

-0 


1  - 1 


(3.6) 
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B\;j  :  tangential  xeiociiy  component  induced  at  the  i''^  panel  control  point  by  unit 
strength  \onicity  distribution  on  the  j^"'  panel 


Bb  -  cos  (8  -  0  )  |3  -  sin  (0  -  0  )  In 

11  *1  ^  !  i  ^ !  1  1 


i|  -  ) 


J.7Z  1  i  U  .  :  p 


(3.7) 


!|  _ 


=  1/2 


'.vhere  the  geometncal  quantities,  depicted  in  Figure  (3,3)  are  defined  bv 


(\ni  -  ~  -  (ym  -  y  ) 

%■  -  x 


\m 


I  I  ^ 


vm. 


I  -  :  ~  } 


0  -  arc  tan 


V  ,  - 

•  I  i 


(3.8) 


-  '^1 


p..  =  arctant  --  -  '  -  arctan  I  — - 

’’  vxm  -X.  ,y  'xm  -xj 

The  first  boundary  condition  requires  flow  tangency  at  the  control  points  ; 

(V"  ),  =  0  .  i  =  I,  2 .  n  .  (3.9) 

In  terms  of  influence  coefficients  (with  \’^  =  1). 


n  n 

Z  ^  . n  .  (3.10) 

.1=1  J  =  1 

The  second  boundary  condition  is  the  Kutta  condition,  which  states  that  the  pressures 
on  the  lower  and  upper  panels  at  the  trailing  edge  must  be  equal  if  the  flow  is  to  leave 
the  trailing  edge  smoothly. 
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I pnfR‘1 

<) 


1  -  I 


'  '^in,  ,  VIII.) 


panel 


limine  (3.3)  KelaHoiisliip  IK'tweeii  (ieoiiieirical  (>ii:iii(i(ii',s 
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I  'sine  nciiunilli's  C(|unti()ti,  ilic  jircssinc  coclUcicnt  can  be  expressed  as  ; 

r  Ic'ta 

r  -  I  -  ^  (;!  I  I) 

r  1  ^  I  , 

pV  ^  ' 

2 

I  he  pressure  equilii)iiuiu  also  implies  criual  velocilics  lor  incompressible  How.  Since 
Ihenoimal  velocities  are  taken  to  be  zero,  the  boundarv  condition  may  now  be  staled  as  ; 

(  V'  l|  -  (  V'  )„  (.V12) 

'aIk’ic  the  negative  siun  is  strictlv  due  to  the  adopted  convetition  of  positive  tangential 
aehK'ities  in  the  direction  (d'  increasing  node  numbering.  Since  the  How  is  positive  to 
the  fiuht,  the  panel  downstieam  of  the  front  stagnation  point  will  have  negative  values 
for  computational  purposes  only  It  is  important  to  note  llial  tiot  all  the  lower  sutface 
panels  have  a  reversed  sigtt,  only  those  downstream  from  the  stagttation  point.  This  is 
especially  signiftcant  for  non- symmetrical  airfoils  or  any  airfoil  at  an  angle  of  attack. 

In  teims  of  inlluetice  coelTicietils,  the  normalized  equation  bccotnes 

II  n  II  n 

IJ*  cos(o,  0.)  V  i  A*  (1.  j  7  V  Ij'.  -  cos  (r/.  0  )  (.V13) 

II  I'  ^liih  ^u  I'  ' 

I  -  I  I  -  I 

I'quations  3,10  atid  3  13  icprcsent  a  linear  algebraic  system  of  (  n  1  I  )  equation  and 
(ml)  unktiowns  I'hc  unknowns  are  the  source  strengths  which  vary  from  panel  to 
panel  (  (p  (j„  )  and  the  votticity  strength  7. 


A  (|. 
M  h 


■/I 

I  -  I 
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Hxpanding  and  rearranging  ilquaticn  i'J  tor  an  exampie  aiiioii  oi'  n  luxies  and 
panels  results  in 

A,'|Lj,  -  A'Acj-,  -  -  B'.A  -11'.'  'in((/-‘',) 

...  -  -  B'l,  ...  B'-!  -  -inta  l.i  !4) 

.q,  -  A-:;  ..  -'/(B-,  ,  ii’A  ,  ...  !]A  -  'int'/ 

rite  equations  now  leadii’.  lend  iheniseives  ti>  solution  in  inairix  :  .irni  ixecastinu 
.vith  simpler  notation,  ihe  lerrns  (eoerficients  ot'  o,  i  mav  be  renaineu  .;nd  the 

sum  of  ail  B'V|  terms  in  parcntltescs  (eocfficients  of  v  )  renamed  a. .  'vlicre  i  1. 

1 .  n  and  i  1  .  2 .  n.  The  terms  on  the  right  side  of  ilte  equation  mav 

be  renamed  b,  .  The  (n-r-l)  equation,  or  in  this  example,  the  74''‘  equation  comes 
from  Equation  .i  13  in  a  similar  manner; 

( A'  I  •  .A-,  |)q,  •  (A‘  .  •  aL  Jq,...  ^  ( A'.  ,,  •Al. 

-7  I  (  1,1  -  B'  1  )  •  (  B‘  i  :  IT  ,  )  •  ,  ■  (  B‘  ,  B'  )  ] 

^  cos  (at)])  -  cos  (a-<}-J3)  (3.1.‘^) 

The  coefficient  q,  may  be  renamed  a-j  j  \I1  of  ilic  i3'  terms  in  the  brackets  together 
Ibrm  the  coefficient  of  y,  now  renamed  a  7174  The  entire  rigid  side  of  the  equation 
constitutes  the  new  term  itsj 
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I'inalK  c\prcssinu  this  svstem  can  he  written  in  ^trncise  matrix  form 


,  I  ; .  ■'  I ,  n  r  I  q ,  h , 

1  : .  n-i  '•I2 

(3.16) 

a  ,  1  .  a  a  h 

:i.  1  _  a.  n  •  I  '■I,,  ' 

a  .  ,  a  . a  ,  •/  b 

Tliis  svstem  is  solved  in  the  program  usimj  a  (iausaian  (-limination  subroLiiine.  With 
ihe  values  oi'  the  q,  and  y  known,  the  velocity  at  each  panel  control  point  may  be 
calculated 

it 

-cos(a-O^)  .  i  -  i,  2,  .  n,  (3.17) 

I  .  I 

The  total  velocity  is  equal  to  the  tangential  velocity  due  to  taking  the  normal  velocity 
to  be  zero 
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!).  ^^ST^:AD^  I’ANF.l.  (  ODE 


I'he  cMcnsioii  ot'  mis  iiietnoo  .m  airtcMi  c\i‘'oi'iL'n(.'inii  ..  :'’nm  ^I'lan'jc  in 
mule  aiiack  can  tie  aecompiisiieu  as  ibilnvvs  ii  ii  ;iiiponam  to  lecaii  I'l.matrs 
ean\  ilau  ■  isuaiization  exDenmems  v.nich  liemimsiraiea  ti'.e  ueneraiioii  it  a 
starting  \onex  .irV  the  traiiing-edge  \'-hcnever  the  angle  m'  attack  ivas  citanged. 
i  lence  continuous  citanues  in  angle  oi  attack  i'rocluce  a.  continuous  Mieudini:  ot 
sonicitv  into  tlie  iraiiinu  wake  This  shedding  can  tie  cxpiaineu  i'v  me  lieiinnoltz 
ortex  itieorein  'Aiiicii  icouires  mat  anv  ciianue  in  me  circulation  aiounu  ’.lie  ainbil 
must  be  matched  bv  the  appearance  of  an  ec]uai  counter  \  ortex  ..t  me  uaiiing 
edge  in  order  to  acitieve  constanev  of  tlie  total  circulation  in  me  :iou  ii.eld. 

The  vortices  shed  from  the  trailing  edge  move  with  the  tluid  panicles  of  the 
surrounding  Iluid  and  hence  are  swept  downstream  with  a  speed  essentially  equal 
to  the  free-stream  speed  Therefore  these  voniccs  will  stav  close  enough  to  the 
airtbil  for  a  finite  time  to  influence  its  pressure  distnbution  T.nv  unsleadv  ainbil 
theory  therefore  must  describe  this  \onex  shedding  process  ii  is  tliis  feature 
whicn  distinguisiies  unsteady  ainbii  iheorv  from  iis  sieadv  counterpart  V!ie 
above  considerations  suggest  that  the  airtbil  can  be  modeled  by  similar  source 
and  \oncx  distributions  as  in  the  stcadv  case. 


‘Phiivcr.  .M 


CIliss  i  cciiiic  Nnic-.s  Mioiil  I’osier.'uiliiiic-  scIpkiI,  .Mnincic'. 


.lllldllllll 


42 


Hence  source  panels  are  auain  placed  on  the  airi'oii  surl'ace  together  a 
sinuie  \onex  distribution  which  is  the  same  ibr  each  panel.  .\n  arbitrary  cnange  in 
angle  of  attack  then  is  subdivided  into  small  step  changes  such  that  a  .staning 
\ortex  is  shed  into  the  wake  each  time  a  step  change  occurs  and  the  tlow  tangency 
and  Kutta  conditions  are  enforced  at  each  time  step. 

In  addition  to  the  n  airfoil  panels  awake  panel  is  assumed  to  be  attached  to 
the  trailing  edge  such  that,  after  each  time  step,  the  sorticitv  of  the  wake  ;:anel  is 
assumed  to  be  concentrated  into  a  single  point  sortc.x  which  detaches  trom  the 
airfoil  and  stans  moving  downstream  with  the  free-stream  selocitv  \s  .'iie  wake 
panel  detaches  from  the  airi'oii  a  new  one  is  created  and  a  wake  ot  point  vortices 
is  created.  As  in  the  steady  problem  ii  unknown  source  strengths  and  one  unknown 
vorticiiy  strength  on  the  airi’oii  are  introduced. 

In  addition,  there  is  now  the  unknown  vorticity  strength  of  the  wake  panel 
whose  length  and  orientation  are  also  unknown.  Hence  there  are  three  additional 
unknowns  which  need  to  be  determined  by  a  proper  set  of  equations;  (a)  the 
How  tangency  conditions  provides  ii  equations,  (b)  the  Kutta  condition  provides  an 


aUUlllUllCll  ^^UallOlL  llsv  i  i;ivvyiviii  uii  wvjuuli\/li 


ity  strength  on  the  wake  panel. 


Therefore  two  more  equations  are  needed  to  determine  the  length  and  orienta¬ 
tion  of  the  wake  panel.  These  are  obtained  bv  making  the  following  assumption.s.  i  e  . 
(d)  the  wake  panel  is  oriented  in  the  direction  of  the  local  resultant  velociiv  at  its 


nid-poiiii.  ;e)  the  iciiuin  the  ■/.ake  panel  prorionionai  te  ihe  niauimucie  la  the 
'.v)eal  resultant  \eicK'itv  at  us  iiiid-poini  and  to  die  iiine-sicp  I'iie  ;\uita  e.u'uiition 
and  the  last  two  ccindiiions  arc  nonlinear  equations  .iiid  iherei'ore  necessitate  an 
Iterative  solution  procedure  i  iuure  i.vdi.  from  ilef  i25j.  shows  the  panei  metnod 
representation  for  unsteady  How 
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I Icimholiz  s  theorem 


Vortex  Slietlditig  nt  I  itiic  Step  t^. 


'V  (AJk  '  f'k  *'k. 


I’nitci  j 

Source  I )i$(ribiitioii 
Voriicity  iJistrihution  y. 


j  i»  I 


k-2  "  *  k-l^ 


^  A. 


^  *  t.  -  J  "  *  k  -  2  ^ 

C^) 

*  *  k  4  "  *  k  •  J 


Iaii0  -  — 

(UJ, 


^  \o. 


\  -  (»k-‘k.i)yimv)k^  ’  (^X'l 


I'iguic  (3.1)  I’niirl  Mcthorls  Ucpi rsriitalioii  (or  (JnstcnrJy 
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rHEODORSKN'S  ANAIASIS 


As  stated  beiore.  i’lieoaorsen's  tiieorv  lor  oseiilatinu  ii.u-piate  airi'oiis  in  iiixiscid 
incomnrcssirle  ''.ow  is  \'-eil  .iocumentea,  '^hereiorc.  ■  e  ratnimanzc  .  ore  ".iv  the 
llnai  formuias.  Fite  nit  on  an  airfoil  osciilatinu  at  tfeaticncv  is  enen  r\ 


t  p  (■)■  ''  ’  ;  i.:.  I  it.,. b)  I  1..,  -  I  !  2  1  )  i.;. 


181 


nero 


Li,  ^  .  2  i  ( '(k)  )  k  . 

L,  -  1/2  -  (  i/k  ).  1  i  -  2  C(k)]  ••  (2  /  k-  )C'(k). 
k  =  '.0  b  /  L'  (  b  is  the  semi-chord  and  "a"  is  the  clastic  axis  i 
C(k)  is  the  lift  detlcicncv  function  and  is  utven  by 

H'*’  (k) 

C(k)  --  . - ^  -  -  =  !■  (k)  -  I  U(k)  (3  19) 

Il|-’  (k)  -iH*-’  (k) 


'.vhere  Ih.dk)  is  the  Hankel  function  of  second  kind  and  is  uiveit  in  terms 
of  I3essei  functions  as.  ll“„(k)  J^lk)  ^^,ik) 
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F.  GARRICK'S  PROPULSIVE  FORCE  ANAL\S1S 


Gairick  15]  applied  Theodorsen's  tlieor\'  to  determine  the  thmst  ,ind  drag 
generated  bv  an  airtbil  which  executes  a  sinusoidal  motion  in  pitch  or  in  plunge. 
He  found  the  follow inu  formula  for  the  average  propulsive  torce 

a,h‘  -  (a,  -  -(a,  -b,)(/  li^  I.V20) 

where  a|  !'“  Ci' 

a.  b-  :  a,  i  I  k-  ■  I  1 '2  -ar  ]  -I  4  -  1 1 ’2  -  a  iF  -  (1  k  : 
a4  -  b  1  a|  [  -  sin((p2-<Po)  -  d'-  -aicos((p2  -  <h,,)  1 
-  (F/2)  cos  (  (P2  -<Po)  ~  sin  ('  (P2  -  <Po)  1 
b2  -  b“  [  -  ay2  -  F/  k"  +  ( 1/2  -  a  )  (G/k)  ] 
b4  (b/2)[  (1/2  ^  G/k  )  cos  (  (p2~  )  -  (F/k)  sin  (  (P2  -  <Po)] 

For  pure  plunge,  the  average  propulsive  force  is  reduced  to 

P^  -  Tipbco^h^  F“(k)4G‘(k)  j  (3.21) 

For  pure  pitch,  the  expression  for  the  average  propulsive  force  is  given  by  ; 

P,^  -  Trpb’«“a“  [a,  •  b.,]  (3  22) 

where  an  and  bi  are  defined  in  Equation  (3.20) 
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The  expression  for  ihe  iuerage  drau  in  1-quaiion  (.'  20)  '.^f  (jarricK's  leoori 
'o7  1 5], was  found  to  he  in  en'or  The  correct  cxnression  as  uwcn  in  tjarrick'^ 
subsequent  .cview  iOO]  is 

-  1  1  "  i  I  1  (i; 

i:,  ,,  V  C“  -V  -  -V  .  ..v23) 

■'  'I  .  ,  .  ..  -  ,  ) 

k"  "  -  -  s'  -  ■  ■ 

A  here 

('"  t'"  -  (j-  .aid  \  -  .T  .  |..  (•)'■.  /,■  ' 

G.  (  OMPARISON  OF  RESULTS 

In  a  recent  thesis.  Riester  (  43  |  presenteu  a  uetaiied  comparison  ■.ith 
Theodorsen's  predictions  '  the  lift  and  pitching  moment  computed  with  the  un¬ 
steady  panel  code.  Generally,  good  agreement  was  found,  thus  lending  credence 
to  the  reliability  of  the  panel  code.  In  the  present  work,  we  are  interested  in 
the  prediction  of  the  thmst  and  drag  forces  generated  bv  an  oscillating  airtbii. 
Calculations  were  performed  to  ascenain  the  panei  code’s  capability  .  Figures 
(3.5)  to  (3.7)  show  the  comparison  between  Garrick's  predictions  and  ilie  oanei 
code  results.  These  figures  present  plots  of  the  average  propulsive  force  coef¬ 
ficient  vs  the  non-dimensional  plunging  amplitude  of  oscillation.  HB.AR  (IfBAR 
=  h/c  ),  for  three  reduced  frequencies.  k  =  U.l  .  1  and  2.  The  two  methods  are 
in  good  agreement  for  small  plunge  amplitudes,  ks  might  be  expected,  the 
agreement  starts  to  deteriorate  with  increasing  plunue  amplitude.  The  panei  code 
was  applied  to  an  \ACA  0009  airfoil  using  100  panels. 
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H.  KS  l  I.MA  i  l  [)  POWI  R  KKI)II( HON  DI  K  TO  lUK  !  OR  on-6A 


HK.LK OPI  I  R 

i’ho  Lillnnaic  (Mijeciix'c  of  tliis  research  'Aas  k'  ojiciiiiine  tiic  eii'cci  oi  iinsteadv 
aerocivnamics  on  tlie  measured  power  \ allies  .ai  tne  OH-OA  lieiicomer  'v’e  are 
especiallv  interesied  m  knowing  '.'.helhcr  these  eii'ecis  couiti  I'roduce  the  power 
■jiianges  recorded  duiinu  ;!ie  NASA  -  Army  imen  l.'on  I’ighl  lesis  in  the  lliuhl 
rcuime  i'rom  hover  ic-  100  knots 

In  I'i'dcr  10  esiiniaie  the  chanue  m  torque  Uue  lo  propulsive  riMce  i  loot- 
pounds  per  Itlauei  due  to  HHC'.  ihe  hlade  span  was  uiviued  into  sections  \s- 
sLiming  liover  conditions,  itie  reduced  frequency  at  eacii  hlade  section  was  cal¬ 
culated  The  panel  cttde  was  used  to  compute  the  propulsive  force  coefficient. 
Cd,  at  each  blade  section  assuming  pure  plunge  motion  with  a  plunging  am¬ 
plitude  of  I",  fhe  blade  deflections,  obtained  from  Wood  et  al  (27].  were  then 
used  to  estimate  the  real  value  of  the  plunging  amplitudes  at  each  blade  segment. 

Finally,  the  net  blade  torque  could  be  calculated  bv  considering  the  propulsive 
i'orcc  at  each  blade  section  and  multiplying  this  vaiue  limes  the  lotor  radius,  men 
summing  over  the  length  of  the  blade.  Table  3.1  .  from  Ref  [42],  shows  the 
()H-6A  helicopter  parameters  that  were  used  lo  estimate  the  piopulsive  Ibices. 
The  most  imporiant  parameters  are.  number  of  blades.  Q  -4.  lotor  radius  (ft). 
R  ^  13.2.  blade  chord  (ft),  c  ~  2b  '  0  57.  rotor  speed  (rad/sec),  il  -  4'4,2.  and 
helicopter  weight  (lb),  (j  '  2550. 
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rABL[:3.l  HIE  ()H-6A  helkopter  rotor  parameters 

Rotor  ladius  (ft)  la  2 

Number  (U'  blades  -I 

Rotor  ■■•peed  ( rad/sec i  'O  26 

Ulade  chord  dti  6  57 
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Table  (3.2)  shows,  at  eacn  blade  section  r.  R.  the  reduced  iVcquency.  k(r). 
local  speed.  UCr),  the  computed  average  propulsive  ibrce.  Py,  torque. 

T.,jjW'i.  due  to  the  propulsive  ibrce  at  each  biade  -ection.  'umming  me  esti- 
atated  torque  at  each  blade  section,  the  resuiiani  tcroue  aue  to  the  pronuisi\  e  rbrce 
per  blade  was  hound  to  be  8,!i  i't-lb  assuming  pure  oiunging  motion  with 
constant  plunging  amplitude  of  1  "  at  each  biade  -  action,  ''.hsod  et  ai  12h’].  ap¬ 
plied  Garrick's  equation  and  considered  the  effect  ca'  the  second  to  twelfth  har¬ 
monic  loads.  He  found,  for  2P  harmonic  (  ITHC  "  n  '  '30  deuree  nnase  at  OO 

knots),  that  the  power  gained  is  about  11.3  horsepower.  Tor  more  detail,  the 
reader  is  referred  to  laole  2  of  Ref  [271  Figures  .'  3  to  i  '  :.now  me  oiade 
detlections  for  several  harmonics  used  bv  iliis  reference  -hown  m  Figure  (3. 12) 
IS  the  average  drag  coefficient.  Cdav.  \ersus  reduceu  frcquencv  as  obtained  with 
the  panel  code  for  a  plunging  amplitude  of  1’’  The  ainbil  used  was  \’.A,C.\  0015 
(OH-6A  airfoil).  The  propulsive  force.  P.^i, ,  is  given  by  . 

=  0.5.  p  .  U“  C  dr.  Cj  (3.24) 

and  the  torque.  Txh,  at  each  blade  section  is  given  bv: 

^'.h  ~  f*\h  •'  •  -v2c) 

where  U  is  the  local  speed  at  each  section ,  Cj  is  the  propulsive  force  coef¬ 
ficient  as  obtained  by  the  unsteady  panel  code  and  r  is  the  radial  distance  of 
each  blade  section.  For  the  OH-6A.  the  angular  freauencv .  Q  ^  30.26  rad/sec. 
therefore  the  power  per  blade  is  found  to  be  0,74  horsepower  as  given  bv  : 

P  =  Txh  .  O  .  (3  26) 
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l  ABU:  3.2  I.OC  Al.  REDI  CED  I  RF.Ql  KNO  .  K^(r).  PROIM'I.SlVi: 
i  OR(  i:.  i\h(i  )  AND  roRQri-:.  i\,, (r). 


r/R  k(r/R)  Kr/R)  |  l\h  (r/R) 


1  - 1 

!  -  1 

1  ft  '  sec  1  ; 

1  ibf ; 

0  1 

1 1  iS(i4 

60  4  1 

6  0()18 

1)  2 

'1  4,1 1 

132,8  ' 

( 1  67^*8 

0  .3 

0  2877 

1 09  2 

I)  9530 

0,4 

0,258 

205  6 

0  108 

O.f 

0,173 

332  0 

0  118 

0  1 

0  144 

398,4 

0  127 

0  7 

0,1233 

404,8 

0,134 

O.S 

(J  1079 

531,2 

0  139 

1  )  ’  ) 

i 

1 1  6)05'! 

507  6 

1,'  144 

I\h(r/R) 

;  I't-lb  ! 

I  1:2 
' '  2(' 

0  4'4 
0  64 
0,80 

I  09 
I  42 
1 ,50 
'  8! 


S  1  i'l-lb 
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Figure  (3.8)  Blade  ryellcctiou  vs  Radial  Station  at  2P 
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I'igurc  (3.9)  Blade  Deflection  vs  Radial  Station 
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l  iemc  (3.11)  Flladc  Dencdioii  vs  Radial  Slalioii  al  5P 
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Reduced  F''reiiuency 


I  AIUJ:  3.3  :  (  (l,„  versus  RF.DI IC.i:!)  MUCQUKNfA  K,,  at  libar-  I  " 
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IV.  ANALYSIS  OF  WAKE 


INTERFERENCE  EFFECTS 


A.  INTRODUCTION 

The  results  of  C'hapter  ill  nIiou  that  propulsive  forces,  ueneraicu  -ulely 
Dv  plunac  oscillations  of  an  airtbil  are  insufficient  to  cxpiain  ilte  novser 
!cducii(rn  obser\'cd  and  measured  on  tlie  Itelicopter  Therefore.  ;n  this 

chapter,  've  examine  whether  the  power  reduction  inmht  be  caiised  h\'  I'avorable 
wake  interference  effects  due  to  the  wakes  shed  from  preceding  blades. 

To  this  end,  the  extension  of  the  unsteady  panel  code  to  the  case  of  two 
airfoils  [25]  was  adopted  and  applied  to  the  study  of  the  wake  interference  elTects 
between  two  helicopter  blades.  Wake  interference  effects  were  studied  in  the 
1950’s  by  Loewy  [4]  by  extending  Theodorsen's  analysis.  Loewy  treats  an 
infinite  number  of  layers  of  shed  vorticity  placed  beneath  the  rotor  at  a  given 
wake  spacing  and  dcvelopes  a  modified  lift  deficiency  'unction  His  interest 
was  in  the  possibility  of  wake-induced  blade  llutter.  which  restricted  his  analysis 
to  the  determination  of  lift  and  pitching  moment  changes  due  to  wakes  shed 
from  preceding  blades.  A  closer  inspection  of  Loewy's  paper  suggests  that  it 
should  be  possible  to  use  the  same  approach  to  examine  drag  and  propulsive 
force  changes  on  tlie  rclcreiicc  blade  due  to  wakes  shed  from  nciglibounng  blades. 


60 


I'unhcr.  it  is  desirable  to  earn’  this  one  step  I'uribcr  and  iincsiipaic  the 
ease  oi'  a  sirmle  wake  shed  from  liie  preceding  blade  which  intett'eres  with  the 
'.vake  shed  tVom  the  rcterence  blade’  Tliis  ease  allows  a  direct  comparison 
with  the  unsteady  two-toil  panel  code. 

Also,  there  are  regimes  of  helicopter  llight  where  cmiy  the  first  iaver  of 
shed  \onicitv  is  of  consequence.  In  the  following  sections  we  tust  desenbe  the 
two-foil  pane!  code  and  ('ouch's  modification  of  the  i  oewv  anaivsis  This  is 
I'ollowed  bv  a  presentation  and  comparison  of  ;he  maior  results. 

B.  TNSTEADY  TWO-FOIL  PANEL  CODE 

The  extension  of  the  single  airfoil  code  (  previously  described  )  to  the  case 
of  two  airfoils  requires  no  new  building  blocks.  However,  the  two-foil  analysis 
requires  the  introduction  of  five  frames  of  reference,  namely  two  moving  local 
frames  of  reference  which  arc  attached  to  the  two  airfoils,  two  frozen  local  frames 
of  reference,  and  the  inertial  frame  of  reference.  Furthermore,  it  requires  the 
satisfaction  of  the  two  KuUa  conditions  which  are  coupled  non-linearlv  The 
solution  requires  an  iterative  procedure  to  compute  the  two  vorticip/  distributions. 


’  A  detailed  devclupinenl  id'  this  llicoiT  i-.  djnvcd  and  |>wsciili;d 
in  a  rccciil  S'I'S  thesis  i'l  Mark  C  uueli.  Itel'  |  IO| 
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!l  also  iioccssiiates  the  creation  of  a  subroutine  wiiicii  iranst'ornis  ail  coor¬ 
dinates  in  either  of  two  local  frames  in' relercnce  to  the  elobai  ifame  I'f  reterence 
Finallv.  ii  requires  the  extension  ol'  the  inllueiiee  coet’ficient  coiieein  to  inciude 
the  effect  of  the  second  airfoil  with  its  own  wake  and  it  reciuires  the  introduction 
of  an  additional  inlluence  coefficient,  namelv  tiiat  c>n  die  w.ike  element  jue  to 
the  'Aake  element  from  the  other  airfoil.  If  eacii  airfoii  d-  modelled  -ith 
panels,  tiien  this  produces  a  IN  x  i  N-.i  )  matrix  '.nich  suDseuuentiv  -oived 
Iw  (jauss  elimination 

Pane  j  25]  was  onlv  able  to  \enfy  ihis  code  "c  .oniDarmu  me  es.imputed 
jiressLire  (.iistnbutions  '..ith  jirexious  '.cork  bv  (.jiesinu  .iencc  'i  was  nec¬ 

essary  to  funher  evaluate  the  code  by  eoinpanng  its  output  with  otner  known 
solutions  before  applying  it  to  the  oscillatorv  blade  interference  problem.  These 
comparisons  are  given  in  section  C. 

The  main  parameters  used  in  this  code  are.  the  airt'oil  type,  number  of 
panels.  N.  the  relative  locations  of  the  two  airfoils,  initial  angle  of  attack,  am¬ 
plitudes  of  osciilation  ( pure  pitch,  pure  plunge  or  iioih).  frequenev  oa’  oscillation 
for  each  airibii.  the  pivot  poinis  (for  pitch  cases),  rise  timet  for  ramp  motions) 
These  input  parameters  arc  in  the  file  "fort  1" 

This  program  produces  an  extremely  large  amount  of  output  to  the  screen.  It 
is  convenient  to  write  the  output  to  a  file  on  a  T'nix  based  machine  during 
program  runs.  The  output  llle  can  gel  very  large  for  long  program  runs  (when 
the  time  historv  is  long) 
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For  this  reason,  the  logical  \'ariable  "output’'  was  added  to  the  input  tile  in  two 
modes  "  true"  or  '  false".  When  this  \ariable  is  set  to  "false",  the  screen  output 
is  not  printed  and  the  required  outputs  are  redirected  to  output  files. 

The  following  list  describes  the  output  fles  and  the  data  thev  contain: 
fon.2  :  This  is  for  user  supplied  airtbil  coordinates,  if  desired  . 
fort  .)  ;  This  fie  contains  the  global  coordinates  of  .\F1  at  each  time  sieo. 
fon.4  This  file  contains  the  global  coordinates  of  .\F2  at  each  time  step. 
fon.7  This  file  contains  the  lift,  moment,  and  drag  coefficients  for  both  ainbils. 
fort.8  ;  This  file  contains  the  pressure  coeff.  for  .ATI  at  each  time  step. 
forT.9  :  This  file  contains  the  pressure  coeff.  for  .\f  2  at  each  time  step  . 
fort.  10;  This  file  contains  the  first  airfoiTs  core  vonex  (wake)  positions, 
fort.  11;  This  file  contains  the  second  airfoil’s  core  vortex  (wake)  positions, 
fort.  12;  This  file  contains  the  angle  of  attack  at  each  time  step, 
fort.  13;  This  file  contains  the  y-axis  translational  motion  at  each  time  step. 

C.  LOEWY’S  ANALYSIS 

Loewy  [4]  developed  a  two  dimensional  model  for  the  unsteady  aerodynamics 
of  the  blades  of  a  hovering  rotor,  including  the  effect  of  the  shed  wakes.  Figure 
(4.1),  from  reference  [40],  shows  his  two-dimensional  model  of  the  helical 
wakes.  The  wake  in  the  case  of  a  multi-bladed  rotor  in  hover  consists  of 
helical  vortex  sheets  below  the  disk,  one  from  each  blade. 
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Assuniinsi  sinusviicial  motion,  the  (iownwasM  o'.or  ;!ie  roicrcncc  airloii.  \  ^  and 
its  \(iniciiv  (either  l.oiind  or  Nhcdl  mav  he  'Aiittcn  icspeeii\eiv  as 


where  the  iiarrea  (inantities  are  complex  I'lte  nuiuced  \eiocn\  at  a  point  on 
the  airtbii  resiiltinu  from  an  element  of  \orticitv  .a'  strenetn  ',,1  at  a  eeneral 
point  in  tl'ic  wake  is  eiven  hv  j4] 

(\  i.)dl 

■  ii.i  ' 

d  N'.  ( \ .  I )  ..  ,  4  2) 

l.T  I  (  \  •  ( n(,)  •  (j )  “  (  h)  ’ 

where  Q  is  the  total  nttmber  of  blades  ;ind  h  is  the  s ertical  distance  between 
succcssice  rows  of  corticitv.  ti  is  the  number  of  revolution  index,  and  p  is  the 
number  of  blades  index  as  shown  in  Figured  I.  The  wake  spacing,  h  is  defined 
as  [4]: 


h 


27:u 

bQQ 


(4  .3) 


where  u  is  the  inflow  'clocity  (  lime-averaued  nomial  to  the  disk) 
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1‘igurc  (4.1)  Acrodynnniic  Unsteady  Model  for  a  Mnlti-Bladed  Rotor. 


The  total  downwash  over  the  airfoil  can  be  expressed  conveniently  by  summine 
the  integrals  involving  the  bound  vorticity  and  the  row  of  vordcity  in  the  plane 
of  the  reference  blade  with  those  for  the  rows  of  vortices  below  the  plane  of  the 
rotor  : 


w(x.t)  = 


I. 


'-I 


( V  0^  '  (»(? ' 

(4.4) 
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The  first  two  terms  on  the  right-hand  side  of  this  integral  eauation  are  those 
that  arise  in  the  classical  unsteady  airtbil  theor\'  of  Tlieoaorsen  i  '  I.  the  third 
and  founh  terms  contribute  the  intluence  of  ail  the  vonicitv  below  me  mane 
of  the  rotor  disc.  Loewy  solved  the  integral  eauation  using  Sohngen's  inversion 
formula  [311.  and  obtained  his  lift  deficiency  function,  C  (k.m.it).  "vnich  is  gi\en 
hv  : 

C  (k.m.)t)  =  F  (k.m.h)  -  i  G  (k.m.h)  .  i4  .>) 

where  k.  is  the  reduced  frcquencv  (defined  bv  tb/L'f  tri  is  uwen  iw  the 
ratio  to/Q.  h  is  the  wake  spacing  defined  before,  and  the  real  ana  imaginarv 
pans  of  Loewy's  lift  deficiency  function.  F  .  G  are  given  bv  [4]  ; 

J.aA-  (Y,  T-pJ,)B 

F' (k.m.h)  =  -  -  . -  -- 

-f  B" 


-G’  (k.  m.  h) 


where 


(Y,  ^pJ,).A^J,aB 
.A“  -  B" 


-A  =  J,a  + Y  -J  p 

1  O 


.  -B  =  Y-jp-J  a 


and 


kh  -kh  ^ 

e  -e  ^  2sin(2:tm) 

a  =  - - - -  .  P  =  -  --  . . .  ; 

kh  —  .  —  .  ~kh  kli  — rvii 

e  -  2cos  (2:tm)  -  e  e  -  2cos  (2;tm)  -  e 
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Here  it  can  be  seen  that  as  the  spacing  between  rows  of  vorticity.  h. 
becomes  infinite,  f’  approaches  F  and  G’  approaciies  G.  /Mso  as  k  lends  to 
infinity.  C’tk.m.h)  approaches  Cfk).  When  h  is  zero,  ail  the  vorticity  lies  in 
the  plane  of  tlie  airfoil,  singularities  result,  and  the  meaning  of  any  obtained 
results,  physically  speaking,  is  not  clear  [4],  Loewy  [4],  and  Couch  [40]  showed 
plots  for  F  and  G  versus  reduced  frequency  k.  Loewy  treats  an  inllnitc 
number  of  wakes  and  Couch  explores  a  finite  number  including  the  special  case 
of  a  single  wake.  I  hcsc  plots  show  the  effects  due  to  wake  phasing,  m  ,  and 
wake  spacing,  h.  A  case  of  special  interest  is  the  wake  sheddiruz  from  one 
jMcccding  blade  only.  For  this  case  l.ocwy's  lift  deficicnev  function.  C'  .  ;cduccs 
to  the  following  lift  function.  C*,  obtained  by  Wood  and  Couch  [40]; 


where 


C‘ 

F* 


=  F*  -  i  G  • 
/|(l 


(4.6) 


An  . (1  - 2i,„)  *  r„  -  21, ^  B,---  j'  -  2.;^ .  ^ 

,v 


I 


The  propulsive  force  coefficient  can  be  calculated  and  the  results,  obtained  by 
Wood  and  Couch  [^0],  using  F*  and  G*  are  ; 


For  pure  plunge 
For  pure  pitch 


Jc 


-T 


1/1  ' 


. 
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1).  i:\  ALrATION  OF  THE  lAVO-FOII.  PANEL  (  ODE 


'^e\ci'ai  anaivses  aae  axaiiable  m  liie  a.cmuviiainic  iaeraunv  .'.nicii  I'crniii  a 
niuea  evaluation  of  ilic  iwo-foii  panel  eoue.  i  hcbc  arc 

"round  Effect  on  Airfoii  Lift  in  Steady  I'lialit 

S  .  ika  el  al  (44|  analv/ed  llie  eround  eiTeci  on  a  I-I)  tlai  niaie  '.r'.inu 
ihc  nioihou  of  eoa.foiTnal  iransformaiions  lonure  ^I  d  sliov.s  me  panel  code 
resLiils,  f  inure  4  Ir  icproduces  i  omoiika  resuiis.  ii  Ncen  mai  me  i.'ao  metnods 
are  m  nood  aureeineni.  lauurc  4  4  a  and  Finure  4  4  b  show  liie  comparison  he- 
iween  the  resuii  obtained  by  \Vegiey  |4o]  for  a  i'.)“o  thick  KALlbl  ain'oil  in 
ground  etVect  and  the  pariei  code  prediction  for  a  NAC/\  oOlO  at  zero  angle  of 
attack  The  two  results  are  again  seen  to  be  in  good  agreement  considering  the 
fact  that  the  RAElOl  airfoil  (whose  coordinates  are  not  known  )  was  approxi¬ 
mated  by  a  NACA  0010  airfoil 

2,  flio . e  Effect 

Glauert  [45]  defined  a  factor  B  which  gives  the  lift  reduction  on  the  airfoil 
in  an  unstaggered  biplane  configuration  compared  to  ihe  lift  on  a  monoplane  at 
the  same  incidence.  Table  4  1  gives  the  correction  factor  for  several  non-dimen¬ 
sional  biplane  spaemgs  It  is  seen  from  this  table  and  from  i'igure(4  4)  that  the 
two  results  are  again  in  reasonable  agreement. 
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IABIJ:  4.1  ( OKKIXTION  lACIOK  I  OK  A  IJIIM^ANI 


_ '  n  (  rodcn  n  ( icr.  |45| )_j 


0.50 

0.770 

0.730 

0.57 

0.847 

0.800  ! 

1. 00  i 

0.902 

0.855  : 

1. 25 

0.942 

0.895 

1. 50 

0.97 

0.920 

Another  gunind  interrcicnce  study  is  due  to  Bagiey  [46],  As  can  be  seen  from 
tlie  comparison  of  the  panel  code  (Figure  4  5a)  with  his  predictions (  Figure  4,5 
li  ),  in  ground  clicct.  the  lift  is  reduced  over  the  upper  surtace  and  increased 
over  the  lower  surface  As  the  airfoil  approaches  the  ground,  the  flow  between 
the  airloil  and  giouiul  is  reducetl  and  hence  the  pressure  is  increased. 
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Cp  vs  X/C  (SS.Wc=.>.  AOA--  I  1  cleg.) 


Mguic  (4.5  ii)  Piossmc  Distrihiilion  vs  X/C  (  Panrl  Code  ) 


I'ifiiiro  (4.5  l)>  |» 


ressiiic  fXsJnfuMion  vs  X/C  (  U.irIcv) 


I'.  \N  \L\SIS  OF  OSCM.I.A  rOin  liLADF.  < M  FKFFRKNC  F  i  FFIA  TS 


iiavinu  Jciermined  in  'v  haptcr  i!l  ihai  pure  ilanpinu  a  inme  '^eiieontcr 
rotor  blade  is  unlikely  to  produce  a  drau  reduction  Miit'ic'cm  to  explain  the 
()H-oA  HHC  iXHver  reductions  [  1  1.  bolelv.  aue  :  ■  me  Kat/.niasr  miect.  .e 
investigate  in  thts  section  the  possibility  ok  lavoraDle  biaue/wake  interi'erence 
etTects.  Consider  the  arrangement  ok  tlte  t\vo  airkoils  'itown  m  i  iuure  i-^o)  phe 
wake  shed  kroin  airkoil  ■■2  impinges  turkoil  ••  I.  ;l'  the  \erticai  spacing  I'etween 
the  two  airtbiis  kails  nelow  a  cenain  s aiue.  depending  on  ilie  aiiiniituue  oi' 
oscillation  ok  airkoil  --2  .  This  case  cannot  he  anaKzeu  wiiii  iltc  unsicauv  p.anel 
code  because  lite  vonices  stan  to  penetrate  the  seeonu  ainbii  or  necause  the 
wakes  shed  from  the  two  airtbiis  come  in  contact  with  each  other,  as  siiown 
in  Figure  (4  7). 

On  the  othei  hand,  if  the  venical  non-dimensional  spacing  is  two  or  more 
no  problems  occur.  A  second  important  parameter  is  the  horizontal  spacing  be¬ 
tween  the  two  airkoils.  This  parameter  controls  the  phasing  between  the  two 
wakes  Hence  the  unsteady  two-lbii  code  permits  the  analysis  ok  I.oewv A  ar- 
rangement  for  the  special  case  ok  two  wakes.  This  case  is  considered  bv  Wood 
and  Couch  [40]  in  closed  form  and  presented  in  Ref.  [.Id].  Furthermore,  the 
single  airkoil  case  is  recovered  if  the  vertical  spacing  is  chosen  to  be  large 
because  the  effect  of  the  lust  airfoil  on  the  second  must  diminish  with  increas¬ 
ing  distance  between  the  two  toils 
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.)/A 


I'ijiiirc  (4.6)  Wake  I’aticrn  (Airfoil  ft  1  IMnngitif'  al  (lish  reduced  I'requi-iicy 
and  (  lose  to  Aiifoil  il  I  ) 
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I•i^ure(4  8).  and  I'iguic  (4,4)  show  the  computed  lift  and  drau(thmst)  co¬ 
efficient  for  airi'oil  !  in  the  presence  of  airfoil  •-  2  Both  airtbiis  are  oscillating 
in  plunge  with  an  amplitude  of  •>  14  and  a  reduced  frequency  of  The 

\errical  spacinti  between  the  two  airfoils  is  200  In  this  case  the  lift  and  thrust 
coefficients  are  found  to  be  identical  to  the  ones  obtained  from  the  singic  airfoil 
code,  thus  showing  that  a  distance  of  200  is  sufficient  to  recover  single  airtbil 
results. 

In  Figures  (4  10)  to  (4  21).  on  the  other  hand,  the  venical  spacing  is  reduced 
to  2  For  these  calculations,  the  time  step  is  .i  to  .f  ’’me  units.  Tlie  plunging 
amplitude  is  again  0  14.  the  reduced  frequency  is  0,0617,  The  twelve  figures 
4.10  to  4  21  show  four  different  cases  of  phasing  between  the  two  wakes.  In 
Figures  (4.10)  to  (4,12),  Tie  wake  pattern,  lift  and  drag  are  shown  for  zero 
phasing,  m  =  0.,  Figures  (4.  13  )  to  (4.15)  show  a  phase  of  m  =  0.208,33,  Figures 
(4.16)  to  (4  18)  show  the  case  of  m  -  0.25,  and  Figures  (4.19)  to  (4  21)  show 
the  case  of  m  =  0  5  , 

.A  closer  inspection  of  these  figures  reveals  the  strong  etfect  of  the  wake 
phasing.  For  m  =  0,  airfoil  If  1  experiences  a  reduction  in  thrust  compared  to  the 
single  foil,  whereas  for  m  =  0,5  airfoil  #  1  experiences  a  significant  increase  in 
thrust.  These  results  are  summarized  in  Figure  (4  22)  where  the  single  airfoil 
results  are  also  shown  as  well  as  the  predictions  by  Wood  and  Couch  [401 
when  specialized  to  the  case  of  two  interfering  wakes  only.  TTie  agreement  be¬ 
tween  the  two  predictions  is  encouraging. 
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xlO-^  Propulsive  Force  Coefncieni  -  Plunge  Only 


0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 


Frequency  Ratio  (m) 


Figure  (4.22)  Average  Propulsive  Force  Due  lo  Plunge  Versus  Wake  Phasing,  in 


Similar  rcsulls  are  ('iblaincd  for  the  jiiich  case.  I’he  airfoils  studied  are  again 
NACA  0007  foils,  pitching  about  the  quarter-chord  point  with  an  amplitude  of 
I  degree  and  a  reduced  frequency  of  0.0617.  lugures  (4.23)  to  (4.32)  display 
the  same  phasing  vaiues  as  in  the  plunge  case.  These  figures  and  the  summa- 
li/.ing  Figure  (4  33)  show  again  that  values  of  m  near  zero  and  one  produce 
increased  drag  compared  to  the  single  airfoil  case,  but  values  of  tn  near  0.5 
produce  thmst.  whereas  the  single  airfoil  produces  drag.  lioth  analyses,  the  anal¬ 
ysis  by  Wood  and  ('ouch  [40]  for  the  special  case  of  two  wakes  only  and  the 
panel  code,  produce  similar  trends  However,  the  .’uantitative  agreement  is  not 
as  good  as  in  the  plunge  case.  The  precise  reason  is  not  sufficiently  understood 
at  this  time  It  ajipcars  that  (he  resolution  of  (ho  suciion  peaks  over  the  leading 
edge  and  therefore  the  computation  of  (he  thrust  is  n.oie  sensitive  in  the  pitch 
case  than  in  the  plunge  case.  Furthermore,  it  must  l.c  s jmembered  that  Loewy's 
analysis  for  an  infinite  number  of  wakes  and  that  by  Wood  and  Couch  [40]  for 
two  wakes  are  based  on  the  (lat-plate  assumption,  whereas  the  panel  code  requires 
a  minimum  blade  tlnckness  in  order  to  obtain  accurate  results.  further  evalu¬ 
ation  of  (he  panel  code  can  be  done  for  the  drag  produced  by  a  single  airfoil 
pitching  about  the  leading  edge  In  (his  case  the  analytical  u'cory  of  Garrick  (5] 
and  the  computational  results  of  Bosch  [15]  are  available  It  can  be  seen  from 
Figure  (4.34)  that  (he  trends  are  again  in  agreemetit,  but  the  '.ode.  results  deviate 
from  (Jarrick’s  and  Bosch's  results  with  increasing  frequency  Tltcse  deviations 
must  be  attributed  to  (he  geometry  differences  (  flat  plate  versu"  NAc'd.  0007  ) 
and  (he  higher  amplitude  used  in  the  panel  code.  Note  (hat  the  pitchit  g  airfoil 
starts  to  develop  thrust  only  at  a  reduced  frequency  gicater  (ban  aboai  1.8. 
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Figure  (4.25)  Ctl  vs  I  *  (  pilch,  in  =  0  ,  h  =2  ) 
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«  7 


Frequency  Hatic  (m) 


I'iuuio  (4.33)  Avouijic  (liiii:  vs  svako  pluisiiif?  (  pitch  only  ) 
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V.  PROPULSIVE  EFFICIENCY 


AND  FLUTTER  ANALYSIS 

A.  INTRODUCTION 

In  the  previous  two  chapters  it  was  shown  that  a  single  airfoil  which 
pertbrms  tiarmonic  plunge  oscillations  generates  a  certain  amount  of  forward 
thmst  ( Katzmayr  effect )  and  that  this  thrust  may  be  enhanced  by  the  oscillaton/ 
wake  sh<’d  from  the  preceding  blade  if  the  frequency  ratio  or  phasing  m  is  in 
the  range  0.2  <  m  >  0.7  approximately  0.5.  It  was  also  shown  that  a  similar 
enhancement  occurs  for  pitching  blades.  However,  pitching  airfoils  develop  thntst 
only  for  relatively  high  values  of  reduced  frequency,  k  >  0.6  by  Garrick 
results  [5]. 

-  urthermoie.  it  is  well  known  that  a  pitching  airfoil  may  develop  negative 
aerodynamic  damping  and  thus  experience  single-degree-of-freedom  flutter.  As 
nown  by  Loewy  [4],  the  possibility  of  flutter  is  enhanced  by  wake  interference. 
This  is  an  effect  known  as  wake-induced  flutter.  It  is  therefore  of  interest  to 
study  the  effect  of  various  parameters  on  propulsive  efficiency  and  on  flutter 
instability  in  more  detail. 
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B.  PROPULSIVE  EFFICIENCY 


The  propulsive  efficiency  is  defined  as  the  ratio  of  the  average  work  of 
propulsion  to  the  work  required  to  maintain  the  oscillation.  If  this  factor  is  zer.? 
or  negative,  we  expect  energy  to  be  extracted  from  the  air  and  fed  to  the  o;  cil  - 
lating  system.  For  a  fiat  plate  performing  a  flapping  motion  given  hv  hit),  where 
h(t)  =  h,,  e''”'  the  result  for  the  horizontal  force  per  unit  span,  averaged  over 
one  cycle,  is  given  bv  [42]  . 


D  =  -Ttpbco'h”  !■“  (k)  +  G~  (k)  (5.F 

where  F(k)  and  G(k)  are  the  real  and  imaginary  parts  cf  the  Theodorsen  lift 
deficiency  function  C(k).  For  the  case  of  pitching  motion,  the  average  horizontal 

force  per  unit  span  length  is  given  by  the  following  equation.  (5.2): 

^  .3  2  2  f|^2  „2'T  1  f\  V'  1  ^  ri  ^  F  Gfl 

D  =  -Ttpb  (0  a  i  F  +  G  -  +  -  a  +  -  F  -  a  --  -  +  a 

k“  u--- 

where  a,,  is  pitch  amplitude  and  (a)  is  the  position  of  axis  of  rotation.  a  =  -  1 
corresponds  to  th';  leading  edge,  a  =  r  j  is  the  trailing  edge.  The  propulsive 
efficiency  then  is  given  by  : 


\V  is  the  average  work  required  to  maintain  the  oscillation.  Foi  pitch,  it  is 
given  by  [42]. 


W  =  7rpb\o'ua;{^(^-a)-(m  J  f(' 


^  ^  r  r  ’  , 
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I'di  pltiimc,  it  is  ui\t'ii  hv 


2  2  t  I 

W  -  Tiph  li  0)  l'(k) 

•’  k 


IIk'ii  the  cITicicncv  factor. 
I''(k)  I  (i'  (k) 


Rpbh' C)"!/!'  (k) 


can  be  expressed  as.  for  plunuc 


(5..M 


'I, 


f(k) 


(.V6) 


and  foi  pitch 


•>  ■  I  '  \  ' '  I  ( 

f  .  (i',  .  i  aj  '  ' 

k‘  “  “ 


!■  (i 
.2  k  ^  2 


(i 


'f  al  -fa.  'Vrf  a,  r 
2''2  ^  ''  2^i  ^2  ■'  kl 


(5,7) 


I'ipiirc  (5  1)  shows  the  inopiilsivc  cnicicncy  of  a  single  plunging  tiiifoilas 
a  (unction  of  reduced  (iec|ucncv  Slowly  flapping  aiiloils  reach  very  high  values 
ol  officionev,  appioachiim  one  As  the  frequency  of  flapping  is  increased,  (he 
elCuicncv  falls  (o  50  percent  hitching  cairfoils,  on  the  other  hand,  start  out  with 
an  efdciencv  of  no  moic  than  ‘^0  percent  at  high  frequency  of  oscillation  which 
iluickly  (alls  lo  zero  and  (o  negative  values  with  decreasing  ficqucncy.  A  negative 
piopulsivc  efficiency  me.iiis  that  no  work  is  required  to  maintain  the  oscillation, 
hut  instead  energy  is  exiiacicd  from  the  airslrcat  Ihis  implies  a  loss  ol  damping 
and  thus  the  possihililv  of  singlc-dcgrec-of-frccdom  flutter  in  pitch 
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Figure (5.2)  shows  the  propulsive  efficiency  of  an  airfoil  which  pitches  about 
the  leading  edge.  The  efficiency  quirkly  drops  to  negative  values  as  the  fre¬ 
quency  is  reduced.  .\s  seen  in  Figure  (5,2),  infinitely  large  negative  values  are 
reached  as  the  inverse  frequency  approaches  24,  corresponding  to  a  reduced  fre¬ 
quency  of  0.038  This  is  the  flutter  condition,  as  will  be  explained  in  more 
detail  in  the  next  section. 

Similar  trends  are  found  for  other  pitch  axis  locations,  as  seen  in  Figure 
5  3  (  pitch  axis  is  at  the  quarter  chord  and  mid-chord  point  ),  The  three  quaaer 
chord  point  is  an  exception.  Here  the  propulsive  efficiency  is  minus  one  at  very 
high  reduced  frequency,  which  falls  to  larger  negative  values  with  decreasing 
frequency  (  Figure  5  4  ).  However,  it  is  interesting  to  note  that  a  small  shift  of 
the  pitch  axis  of  only  0.5  %  in  either  airection  restores  the  efficiency  to  0  5 
for  large  frequencies,  from  which  it  starts  to  drop  with  decreasing  frequency. 

These  results  show  that  there  is  a  fundamental  difference  between  airfoil 
plunge  and  pitch  motions.  Plunge  oscillations  of  a  single  airfoil  always  produce 
thrust,  and  the  efficiency  decreases  as  the  frequency  is  increased.  Pitch  oscil¬ 
lations,  on  the  other  hand,  produce  thrust  only  at  relatively  high  values  of  re¬ 
duced  frequency,  k  >  0.6  This  possibility  of  single-degrce-of-freedom  in  pitch 
flutter  is  discussed  in  more  detail  in  the  next  section. 


110 


12  3  4  5 


1/k 


Kigiirc  (r..3)  rro|Miisive  Efnciciicy  n  versus  1/k  for  an  airfoil  pitcliing 
aboui  Ike  <niar(cr  -  cliord  point  and  nnd- chord  point 
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0.5 


Fipurc  (5.4)  rropiilsivc  ICfficiciicy  )i  versus  l/k  for  an  airfoil 
piU'hiii^  about  the  tlircc  -  (piartcr- chord  point 


C.  SINGLE-DEGREE-OF-FREEDOiM  AIRFOIL  FLUTTER 


Consider  an  airtbil  hinged  at  us  leading  edge  but  elastically  restrained  from 
rotating  about  this  axis  by  a  torsion  spring  with  constant  K,^  (ft-lb  /rad),  The 
airfoil  is  placed  in  a  low-speed  airstream  so  that  the  unstrained  position  of  the 
spring  corresponds  to  zero  angle  of  attack  a.  Figure (5.5),  The  equation  of 
motion  for  this  single-degree-of-frecdom  system  is  given  by  ; 


I  a  +  K  (/  =  VI  . 

\ 


(5.8) 


Mv  is  the  aerodynamic  moment  due  to  a(t)  and  1,^  is  the  moment  of  inertia 
about  the  leading  edge,  .\ssuming  that  a(t)=  ct„e''”‘-  then  equation  (5.8)  can 
be  written  as 


I 

a 

repb"^ 


-  0. 


(5.9) 


where  is  the  natural  frequency  of  torsional  vibration  and  's  given  by 


(j) 


c  1 


(5.10) 


and  m^,  is  shorthand  for  the  dimensionless  aerodynamic  coefficient  and  can  be 
written  as  [41]; 

M 


m  = 

V 


Tcpb 


^  —  -=M  -  iL  H-  j  +a-t-L, 

4  2  10)1  V  n/  \n  J  li  V? 

(0  a  e  -  - 


(5.11) 
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where  "  a  "  is  the  elastic  axis  position  and  at  the  leading  edge  a  ~  -  1.  ii\. 
is  a  complex  number  and  it  is  alunction  only  of  icduccd  ficqucncy,  k  “  (!)  b/lj. 


f'.quation  (5/))  can  be  split  into  ical  and  imaginai'V  parts. 


Re { m  }  -  Re 


TCpb 


(!) 


t  '  (I)  / 


I 


(5.12  a  ) 


Im  I  m^.  ;  ‘  0  (5, 12  b) 


I'lutlcr  occurs  at  that  value  of  the  reduced  freqiiencv  where  the  imaginan.’ 

I  out-of-phase)  part  of  the  aeiodvnamic  moment  becomes  zero,  provided  that  the 
corrcs|ionding  real  (in-phase)  part  yields  a  non-imaginaiv'  llutter  frequency. 


rigurc  (5.5)  An  airfoil  restrniued  lo  rotaJe  about  its  leading  edge 
in  (wo-diincnsional  flow. 
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Figure  5.6  shows  the  variation  of  the  real  and  imaginary  (  damping )  parts  of 
the  pitching  moment  for  an  airfoil  pitching  about  the  leading  edge.  It  is  seen 
that  the  pitch  damping  goes  through  zero  at  a  reduced  frequency  of  0.038  This 
value  agrees  with  the  previously  found  value  using  propulsive  efficiency  ( energy 
extraction  )  considerations. 

Figure  5.7  displays  the  comparison  between  the  pitch  damping  coefficients 
(note  that  the  pitch  damping  coefficient  is  defined  by  Theodorsen  [3]  as  : 

C  =  M  -i'-ai(L  -M,  )--  '-aiL,  .  (5.13) 

The  difference  in  Cm  shown  in  Figure  4.8  is  due  to  a  ditference  between  the 
,  L,, ,  Mj,  and  terms  in  Equation  5,]  3  as  determined  by  panel  code 
(  Riester  [43]  )  and  Tlieodorsen’s  theory  [3],  The  imaginary  part  in  Figure  5.7 
is  obtained  by  Equation  5.11  computed  by  the  panel  code  and  Theodorsen’s 
analysis.  It  is  seen  that  the  two  computations  are  in  reasonable  agreement, 
recognizing  the  difference  in  geometry  (  NACA  0007  versus  flat  plate ). 

Systematic  variation  of  the  pitch  axis  location  shows  that  zero  pitch  damping 
is  possible  for  axis  locations  upstream  of  the  quarter -chord  point.  No  flutter  is 
possible,  on  the  other  hand,  for  axis  locations  at  or  downstream  of  the  quarter 
chord  point  Figure  5  8.  Loewy  [4]  has  shown  that  wake  interference  greatly 
increases  the  possibility  of  flutter.  The  unsteady  panel  code  provides  an  oppor¬ 
tunity  to  investigate  and  to  extend  this  well-known  finding.  This  is  a  problem 
which  should  be  studied  in  more  detail. 
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Real  Gi  Imag.  pirt  of  m  vs  k  at  x  = 


Reduced  frequency  k 


I'imuc  (?.6)  N'iiriiUioii  >vi(h  icdutcd  frcquciiry  K  of  Uio  rr;d  jind  iinnKiiinry 
piiKs  of  rlu'  ooi:-dimt  i,sion;d  acrodynjimic  inoiucid  due  (o  pikliiu«  of 
an  ail  fod  ahoiil  i(s  leadiiiu  <  dt;c  in  an  iiuoinpi  essible  How. 


t  \  0  ij  u  c  G*  u  [re  q  u  e  n  c  y 


I, 

r\ 


I'iuuio  (S.8)  lin:i}>:.  (my  )  vs  U  (  for  ;ui  :iirf<»il  nifrliiiiu  r!i~- 


.5,(>.0,().5,;uul  1.0) 


VL  RESULTS  AND  CONCLUSIONS 


A  DISCUSSION  OF  RESULTS 

This  research  was  undertaken,  to  study  the  etTeci  of  unsteady  aerodynamics 
on  reductions  of  main  rotor  shaft  torque  and  engine  power  measured  when 
the  HHC  system  was  applied  to  the  OH-6A  rotor.  In  Chapter  III  we  study  pro¬ 
pulsive  force  resulting  from  pure  plunge  oscillation  of  an  airfoil  in  invi.scid, 
incompressible  How.  The  results  of  flat  plate  theory,  Theodorsen  [3]  and  Garrick 
[5],  are  used  to  validate  panel  code  results.  Comparison  with  Garrick  showed 
good  agreement  in  pure  plunge,  especially,  for  lower  values  of  reduced  frequency 
and  amplitudes  of  oscillation. 

Figures  3.5  to  3  ,  show  the  average  drag  coefficient  versus  plunging 
amplitude  for  three  values  of  reduced  frequency,  k“  0  1.  1,  and  2  respectively. 

It  is  seen  from  f  igure  3.5,  (  reduced  frequency,  k  is  0.1  )  that  agreement 
between  the  panel  code  and  Garrick  [5]  is  very  good  for  a  range  of  non- 
dimensional  plunging  amplitude  up  to  15  to  20  %  of  the  blade  chord.  On  the 
other  hand,  for  higher  values  of  k,  the  two  results  arc  in  good  agreement  for 
a  range  of  plunging  amplitude  up  to  approximately  10%.  Single  airfoil  code,  was 
used  to  compute  the  time  history  of  the  drag  (thrust)  coefficient  of  the  NACA 
0015  airfoil  (the  blade  section  of  the  OH-6A  ), 
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Figure  3.12  and  Table 3. 3  show  the  average  drag  coefficient  \ersus  reduced 
frequency  for  an  NACA  0015  airfoil  which  executes  plunging  oscillation  with 

1  "  plunging  amplitude.  The  OH-OA  data  presented  in  Table  3.1.  Wood  [42], 
and  Figures  3.8  to  Figure  3.11  from  Ref.  [27],  were  used  to  estimate  the 
propulsive  force  per  blade  obtained  due  to  plunging  oscillation  for  1  P  frequency 
case.  In  these  figures,  the  difference  between  the  deflections  when  the  HHC 
is  "  on  "  and  "  off  "  represents  the  jtlunging  amplitude  The  controller  phase  is 
also  significant  Phase  angle  of  00 or  120”  are  those  where  helicopter  vibra¬ 
tions  were  most  sever,  whereas  at  phase  angle  of  300  ”  or  330”.  the  helicopter 
vibrations  were  significantly  reduced. 

We  look  for  the  greatest  power  benefit  to  occur  at  the  phase  angle  of 
lowest  vibration,  that  is  300”  or  330”.  Table  3.2  shows  the  amount  of  propulsive 
force  and  torque  generated  at  each  blade  segment.  The  results  of  this  table 
showed  that  the  contribution  ol  IP  (first  harmonic)  were  3  horsepower.  Wood 
et  al  [27],  applied  the  Garrick  equation  and  considered  the  effect  of  blade  plunge 
response  in  relation  to  the  second  to  twelfth  haimonic  loads.  He  Ibund  for  2P 
harmonic  only  (  HFIC  "on  ",  33U”  phase  at  OU  knots  ),  the  power  gained  was 
about  1 1.3  horse  power.  For  more  detail,  the  reader  is  encouraged  to  study  Table 

2  of  Ref.  [27]  ‘  . 


'  Wood  ]•  R  Ilit'.niaii.  ,!  ami  Itamcsli  Kolar.  "  I  liulici  I  laniiniiiu  Cmitri'l 
I’roniisc.s  l;np'()\'i.-d  Ih  nainic  Interlace  ( Ipcration.s'  A(iARI)  I'lncecdiniis  <>l 
llic  7S  til  1  liphi  Mechanics  I'anci.  May  |vV|,  Seville.  S|)ain 
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Chapter  IV,  considers  the  influence  of  the  rotor  wake  interaction  and  examines 
whereas  such  interaction  might  cause  favorable  effects.  First,  the  case  of  an  iso¬ 
lated  airfoil  (  single  airfoil  oscillating  in  plunge)  with  absence  of  wake  vvas 
considered  ( or  thr  wake  located  at  great  distance  from  rotor  ).  This  is  the  "Katz- 
mayr  effect"  case  where  the  resulting  drag  force  is  a  propulsive  force,  It  was 
the  primary  case  considered  by  Garrick  [5] 

The  advantage  of  the  Panel  Code  is  that  it  provides  us  with  time  histories 
of  unsteady  lift,  drag,  and  the  wake  trajectory  by  which  we  can  better  understand 
and  analyze  the  problem.  In  each  case,  the  wake  pattern  was  plotted  in  order 
to  be  sure  that  the  wake  produced  by  the  second  airfoil  would  not  inteifere  with 
the  reference  airfoil.  Also  the  lift  and  drag  (or  the  propulsive  force)  were  plotted 
versus  non-dimensional  time,  T*  (  T”"  =  t.U/c ). 

The  closed  form  solutions  are  based  upon  a  very  thin  flat  plate  airfoil 
undergoing  very  small  motions  when  we  model  it  in  the  unsteady  panel  code. 
For  the  calculations  an  NACA0007  airfoil  was  selected.  The  aitfoil  sets  at  zero 
mean  angle  of  attack  and  the  reduced  frequency,  k,  is  =  0,0617 

For  the  OH-6A,  the  value  of  reduced  frequency  ,  k  at  70  %  radial  station  , 
was  calculated  and  found  to  be  .1234  for  the  32  Hz  exciting  frequency.  The 
reference  length  used  by  the  panel  code  is  the  blade  chord  or  2  b,  so  the  0  1234 
reduced  frequency  corresponds  to  0  0617  per  the  panel  code  definition.  Non- 
dimensional  wake  spacing,  for  single  airfoil  analysis,  was  h  =  200  (  mathemat¬ 
ically  equivalent  to  infinite  spacing,  or  we  say  that  the  interaction  between  the 
two  airfoil  is  very  small)  . 
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The  amplitude  ot'  plunge  oscillation  was  It,,  ^  <'.14  (  by  the  panel  code 
.07,  for  the  same  reason  ).  The  panel  code  results  are  compared  to  other  refer¬ 
ences  in  order  to  know  the  code  limits.  The  comparison  of  the  code  results  to 
that  obtained  by  Ref.  (44],  is  shown  in  Figure  (4.2).  This  figure  shows  the 
relative  lift  coefficient  of  an  airfoil  (  a  wing  of  relatively  long  span  )  as  it  ap¬ 
proaches  the  ground  for  different  values  of  angles  of  attack.  The  figure,  showed 
that  ihe  results  estimated  bv  the  panel  code  and  the  results  obtained  by  Ref 
[44],  ate  very  close.  For  h/c  ■  0.5  the  results  obtained  by  the  panel  code 
are  suspected. 

Another  comparison  is  shown  in  Figure  (4.4).  This  figure  and  Table  4.1 
show  the  correction  factor,  B,  versus  the  chord-gap  ratio  in  a  biplane  configura¬ 
tion.  The  correction  factor,  as  defined  by  Ref  [45],  is  the  reduction  of  the  lift 
coefficient  of  an  unstaggered  biplane  compared  with  that  of  monoplane  at  the 
same  angle  of  incidence.  On  Figure  (4.3)  and  Figure (4.5),  the  local  speed  ratio 
and  pressure  coetTicient  vs  X/C  are  plotted  and  showed  very  good  agreement 
with  Wegiey  [46]  and  Bagiev  (1959),  respectively 

The  panel  code  is  suspect  whenever  the  wake  generated  by  an  airfoil 
impinges  on  another  airfoil.  The  numerical  panel  code  method  also  encounters 
difficulty  with  discontinuities.  For  these  reasons  it  is  important  to  carefully  check 
the  time  history  plots  after  each  computer  run. 


123 


In  Figure  (4  8),  tlie  time  history  ot'  the  lift  coetTicient,  C|  variation  is 
plotted.  Shown  in  Figure  (4.9  is  the  drag  coefficient.  Cj  \ersus  the  non- 
dimensional  time  T*.  We  observ'e  from  Figure  (4  8)  that  the  lift  varies  as  the 
frequency  of  oscillation,  with  a  mean  lift  value  of  zero  as  we  expect.  Amplitude 
of  oscillatory  lift  was  found  to  be±  0.0510.  Figure  (4.9)  clearly  shows  the  " 
Katzmayr  effect".  We  noted  that  the  calculated  drag  varies  at  twice  the  plunge 
oscillation  frequency  of  the  airfoil  .Also,  observ'e  that  the  airfoil  is  generating  a 
propulsive  force  as  indicated  by  the  mean  value  of  drag  coefficient,  which  is 
seen  to  have  a  negative  value. 

Looking  to  the  computer  printout  at,  at  the  non-dimensional  time  T*  = 
0.0,  steady  state,  the  drag  generated  by  the  first  airfoil  was  found  to  be 
Cjsjj  =  000131  for  NACA0007.  which  is  the  airfoil  type  that  was  used  for 
most  of  our  analyses  at  zero  angle  of  attack ,  and  a  tO'.al  number  of  panels 
N  =  100  (  i  e50  panels  at  the  upper  surface  and  50  panels  at  the  lower  surface). 
For  NACA0009,  (this  aiP  oil  was  used  to  compare  the  code  results  with  that  of 
Bosch’s  for  a  pitching  case),  the  steady  state  drag  was  found  to  be  .000196 
at  .N  =  100  and  .000113  at  number  of  panels  N  =150 

It  is  known  that,  for  an  airfoil  at  zero  angle  of  attack,  in  an  inviscid 
incompressible  flow  Oeld,  steady  state  drag  is  zero,  so  in  our  calculations  that 
steady  state  drag  was  subtracted  from  the  drag  generated  at  each  time  step.  Shown 
in  Figure  (6.1)  is  the  steady  state  drag  coefficient  vs  angle  of  attack  [37].  Figure 
(6.1a)  shows  the  variation  of  steady  state  drag  coefficient  vs  angle  of  attack  for 
three  different  airfoils;  NACA0003,  NACA0007,  and  NACA0012  at  total  number 
of  panel  of  200. 
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In  this  graph,  it  is  shown  that,  for  the  thin  aitlbil  NACA0003  at  zero  angle 
of  attack,  the  drag  coefficient  is  \ery  small  but  as  the  angle  of  attack  increases, 
the  drag  coeffiicient  highly  increases  more  than  in  the  cases  of  XACAOOO?  and 
NACA0012.  Figure  (0.  lb)  shows  the  effiect  of  increasing  the  total  number  of  pan¬ 
els  N  on  the  drag  coefficient  of  the  airfoil  NACA0012.  As  the  total  number 
of  panels  increases,  the  steady  state  drag  coefficient  accumulated  erroi  decreases. 

Figures  (4  10)  through  Figure  (4  21),  show  the  same  case  as  above,  that 
ot  an  airfoil  oscillating  in  plunge,  but  with  the  introduction  of  a  wake  layer 
in  the  near  proximity  of  the  reference  airfoil  (non-dimensional  w'akc  spacing,  h 
=2  )  for  wake  phasing,  m  =  0,  0  20833,  0.25.  and  m  0.5.  These  correspond 
to  phase  angles  of  0‘*  .  75'*,  OO^’  and  180",  respectively 

We  see  that  the  phase  angles  relate  to  the  phase  relationship  between  the 
reference  airfoil  and  the  wake  layer  of  shed  vonicity.  These  phase  angles  were 
obtained  in  the  panel  code  by  establishing  the  proper  time  phase  between  the 
initiation  of  the  wake  vorticity  and  the  oscillation  of  the  reference  airfoil.  This 
time  phase  is  in  turn  a  function  of  the  frequency  of  oscillation  which  is  directly 
related  reduced  frequency,  k.  For  the  results  presented  here,  for  k-  0617,  this 
corresponds  to  50.92  time  units  for  m-0.0  (O'*  phase);  61.53  time  units  for  m 
=  0,2083  (  75"  phase);  and  76.38  time  units  for  m  =  0.50{  180"  phase),  etc. 
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the  reduced  frequency,  k,  changes,  this  corresponds  to  a  change  in 
frequency  (at  given  chord  length  and  free  stream  velocity  ).  for  k=0  1234,  this 
correspond  50.92  time  units  for  m  =  0.0,  56.23  time  units  for  m=  0.2083,  and 
63.65  time  units  for  m  =  0.5,  etc. 

Figure  (4.10)  shows  the  wake  position  for  in  =  0,0,  the  case  where  the 
wake  shed  from  the  second  airfoil  at  h  =  2  is  in  phase  with  wake  shed  from 
the  reference  airfoil  itself.  In  Figure  (4,21  ),  the  lime  history  of  lift  \ariation 
is  shown 

This  is  indicated  by  the  solid  line  for  ('1|  .  Fhe  dashed  line,  Ch,  indicates 
the  variation  of  lift  on  a  second  airfoil  located  at  the  point  of  initiation  of  the 
wake  layer  at  li  =  2.  Our  discussion  will  focus  airfoil)  ,  indicated  by  the  solid 
line  in  ♦he  figures  mentioned  above.  For  an  airfoil  oscillating  about  zero  angle 
of  attack,  the  mean  lift  is  zero.  We  found  that,  while  the  mean  lift  is  zero,  the 
amplitude  of  lift  oscillation  is  smaller  in  the  presence  of  the  wake  at  zero 
phasing.  From  these  figures,  we  conclude  that  the  effect  of  wake  vonicity  at 
zero  phasing  is  to  reduce  the  "  Katz  nayr  effect". 

Comparison  of  oscillatory  lift  values  of  the  75'*  wake  phasing  case,  with 
the  case  of  no  wake  or  single  airfoil,  shows  that  they  are  identical,  since  the 
values  at  75  degrees  phase  are  the  same  as  that  for  infinite  wake.  Table  6.1 
proves  this  result. 
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The  time  history  of  drag  is  seen  to  vary  as  twice  the  plunge  oscillation 
frequency  of  the  airtbil.  We  also  observe  that  the  airfoil  is  generating  propulsive 
force  as  indicated  by  the  mean  value  of  drag,  which  is  seen  to  have  a  negative 
value.  Results  of  75  degree  wake  phasing  case,  confirm  that  there  is  a  phase 
angle,  m,  at  which  the  effect  of  the  wake  vorticitv  on  the  refeience  airfoil  gives 
the  same  results  as  when  the  wake  is  removed  .  It  has  been  shown  at  a  phase 
angle  of  zero  degrees  that  the  shed  wake  vorticitv'  increases  the  drag  on  the 
oscillating  refeience  airfoil,  i  e.  diminishes  the  "Katzmayr  effect" 

Conversely  we  found  that  at  a  phase  angle  of  about  180  degrees,  there  is  a 
significant  decrease  in  drag  at  the  reference  airfoil  due  to  the  wake  vorticity. 
This  implies  an  intermediate  phase  angle  at  which  the  drag  is  identical  to  the 
"Katzmayr  effect".  It  turns  out  that  such  a  phase  angle  exists.  In  fact,  at  this 
phase  angle  the  lift,  drag  and  pitching  moment  are  found  to  be  identical  to  the 
case  of  no  wake,  and  match  the  "no  wake"  values  identically. 

Table  (6.1)  show  the  aerodynamic  coefficients  (.1(1).  CU2),  Cm(l),  Cm(2), 
Cd(l),  and  Cd(2)  vs.  Non-dimensional  time.  Tfiis  table  is  a  part  of  one  of  the 
output  files.  It  shows  these  aerodynamic  coefficient  vs  time  from  the  nondimen- 
sional  time  unit  113.3  until  159.1.  The  phase  angle  at  which  this  occurs  is  75 
degree  or  m  =  0.20833,  for  the  case  considered  here. 
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liihlc  (6.1)  Acnxiynaniir  coefficients  vs  Non-diincnsioiial  time  (from  time 
step  113  to  159)  for  (he  case  of  wake  phasing,  m  =75  degree,  wake  spacing 
li  =  2,  reduced  frc(|iiency  k  =0.0617  and  plunging  amplitude  hl)ar  =  0.14) 
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0.007188 
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The  proper  phasing  is  achieved  by  positioning  the  reference  airfoil  76.38 
time  units  to  the  rig.  of  the  second  airfoil,  the  airfoil  that  generates  the  wakt. 
We  found  in  the  case  of  the  plunging  airfoil  that  the  wake  at  zero  degree  phase 
increases  the  drag  force,  while  the  wake  at  75  degree  phase  gives  us  close 
results  to  "no  wake  at  all"  case.  So,  we  expect  that  the  wake  at  180  degree  to 
increase  the  propulsive  force.  Shown  in  Figure  (4.22),  is  the  \ariation  of  the 
average  propulsive  force  coefficient  versus  wake  phasing,  m.  The  circles  in  the 
graph  represent  panel  code  results  for  the  parameters  mentioned  above. 

Superimposed  on  the  plot  of  the  figure  are  the  analytical  closed  form 
results  of  the  Loewy  theory  as  modified  for  one  wake,  as  shown  in  Figure 
(4,22)  from  Wood  et  al  [38].  Note  the  good  agreement  between  the  two  sets 
of  results.  It  is  likely  that  the  agreement  can  be  further  improved  by  increasing 
the  number  of  panel  from  100,  used  for  these  calculations,  to  200  or  200  more. 

As  was  done  previously  for  an  airfoil  oscillating  in  plunge,  we  will  now 
consider  the  case  of  an  airfoil  oscillating  in  pitch  with  panicuiar  emphasis  on 
drag  or  propulsive  force.  We  will  first  look  at  the  case  without  presence  of 
layers  of  adjacent  wake,  then  look  at  the  wake  effect  including  the  impoaant 
consideration  of  phasing.  We  will  begin  by  considering  the  numerical  results  and 
comparing  these  results  with  Garrick,  Bosch,  for  "no  wake"  case,  and  then  com¬ 
pare  with  Loewy’ s  theory  when  the  wakes  are  included. 
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Considered  is  the  case  of  an  NACA0007  atrtbil  set  at  zero  degree  angle 
of  attack.  The  airfoil  is  oscillating  with  an  amplitude  I  degree.  The  reduced 
frequency  is  taken  at  k  =  0.06i7  and  the  non-dimensional  wake  spacing  was 
li -200.  which  is  numerically  equivalent  to  having  the  vvake  at  infinitv.  The  time 
history  of  the  lift  variation,  Ci.  and  the  time  historv  of  drag  variation,  Cd,  reveal 
the  following.  We  observe  that  the  lift  varies  as  the  frequency  of  oscillation, 
with  a  lift  mean  \alue  of  zero  as  we  expect. 

In  contrast  to  the  lift  time  history,  the  drag  time  history  indicates  ihat  the 
calculated  drag  varies  at  twice  the  frequency  of  the  pitch  oscillation.  The  am¬ 
plitude  is  given  by  Cd=  0. 000192.  The  mean  cirag  '-alue  is  found  to  be  Cd"= 
0.00009  and  the  mean  lift  is  of  course  0.0  at  zero  angle  of  attack. 

For  an  airfoil  oscillating  in  pitch  with  the  wake  layer  near  the  reference 
airfoil  at  non-dimensional  wake  spacing.  h  =  2,  The  case  of  zero  degree  phase, 
where  the  wake  below  the  airfoil  at  h  =  2  is  in  phase  the  wake  shed  from  the 
reference  airfoil  itself.  Review  of  the  results  indicates  a  reduction  in  the  ampli¬ 
tude  of  lift  variation  for  the  m=0.0  case  when  compared  to  the  case  with  no 
wake  at  all  (wake  at  infinity).  We  see  that  the  lift  amplitude  decreases.  We  also 
see  an  increase  in  mean  steady  drag  accompanies  the  reduction  in  lift  amplitude, 
for  the  case  of  zero  wake  spacing  (m  =  0).  Proper  phasing  is  achic'’ed  at  m 
=0.5.  At  this  value,  w'e  found  that  the  lift  enhancement  in  this  case  of  pure 
pitch,  is  similar  to  what  we  found  in  the  plunge  case. 
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B  CONCLUSIONS 


The  most  important  conclusions  iliai  can  be  drawn  from  this  research  are; 

1.  For  drag,  in  both  cases  of  pure  plunge  and  pure  pitch,  wake  phasing  tan 
increase  or  decrease  the  steady  component  of  drag  acting  on  the  airfoil.  The 
increase  in  the  steady  component  of  propulsive  force  due  to  wake  phasing  in  the 
case  of  pure  plunge  was  sufficient  to  significantly  enhance  the  "  Katzmayr  " 
effect  value  of  [iropulsive  force. 

2.  For  an  airfoil  oscillating  in  plunge  or  pitch,  for  optimum  reduction  in  drag 
and  enhancement  in  vibratory  lift,  the  phase  angle  of  the  wake  vorticity  in  the 
single  wake  case  with  respect  to  motion  of  the  reference  airfoil  should  be  about 
ISO'-’ 


3.  For  an  airfoil  oscillating  in  plunge  or  pitch,  the  largest  value  of  steady  drag 
and  smallest  value  of  oscillatory  lift  occurs  '  'hen  the  phase  angle  of'  the  wake 
vorticity  with  respect  to  motion  of  the  reference  airfoil  is  about  <) '* 

4.  For  lift,  in  both  plunge  and  pitch  in  the  single  wake  case,  wake  pjhasing  can 
increase  or  reduce  the  oscillatory  lift  acting  on  the  airfoi'  Where  in  the  absence 
of  wake  vorticity,  the  effect  of  oscillations  decreases  the  lif  (lift  deficiency  ). 
With  the  wake  present  we  also  observe  lift  enhancement  (  lift  efficiency  ). 

5.  There  exists  a  phase  angle,  75'’  for  the  case,  (reduced  k  =  .0617,  wake 
spacing,  h  2  ),  where  the  effect  of  wake  vorticity  on  the  reference  airfoil  is 
identical  to  the  case  of  the  wake  at  infinity  (  represented  iii  the  code  by  wake 
h  ~  200  ),  or  no  wake  vorticity  at  all. 
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Applying  these  findings  to  the  ()H-6A  case,  vve  can  add  the  following 
conclusions  : 

6.  Measured  reductions  in  power  are  feasible  when  the  "Katzmayr ''  effect  and 
thesis  -  rcponed  wake  enhancements  are  included. 

7.  Measured  higher  harmonic  control  results  show  least  benefit  nt  100  knots 
due  to  the  wake  being  transported  greatest  distance  from  rotor  at  this  speed. 

8  Poor  rcpeatabiliiy  of  measured  open  loop  performance  data  is  to  be  expected 
due  to  shifts  in  vake  position  *. 

9.  For  helicopters,  it  is  standard  rotor  blade  design  to  locate  the  center  of 
gravity,  aerodynamic  center,  and  elastic  axis  at  the  quaner  chord  point  to  avoid 
llutter.  This  was  illustrated  in  the  analysis.  Ilesulis  of  this  thesis  showed  that 
flutter  is  impossible  for  this  condition  without  wake  interference.  However, 
Loewy’s  results  introduce  the  possibility  of  wake  induced  flutter.  Future  investi¬ 
gations.  therefore,  will  have  to  explore  this  possibility  in  more  detail. 


'  li.  It.  WouiJ  Max  1'.  I'liil/.ci.  Ahmed/  .luialiiiia.  Maik  A  (.'uucli 
On  tlic  Onstciidv  Aciodviiainics  ol  llighcr  llannonic  Cniiirol  " 
f'apcr  No.  L'  17,  Ninciccnlh  Im'opciiri  Koluicral'l  I'ORim.  (.  crnolihio 
(C.'oiiioi  llalv,  'icptcinhci  IVV.t, 
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C  RECOMMENDATIONS 


1.  A  new  night  test  program,  directed  at  studying  helicopter  pertormance.  is 
needed.  The  NASA-Army  OH-6A  program  was  dedicated  to  vibration  testing. 

2.  Continued  full  scale  MHC  tests  in  the  NASA  Ames  40 '  x  80  '  wind  tunnel  is 
encouraged.  In  addition,  it  is  suggested  that  tests  be  initiated  in  small  university 
wind  tunnels  to  '  erify  the  findings  of  this  thesis. 

3.  ('ontrolled  whirl  tower  testing  is  also  recommended,  bince  this  allows  careful 
measurement  and  control  of  performance  parameters. 

4.  A  wind  tunnel  program  is  needed  to  explore  similar  cases  and  to  verify  the 
present  research  findings. 

5.  The  panel  method  is  a  significant  analysis  tool  that  can  be  applied  to  study  and 
solve  similar  problems.  It  provides  an  excellent  method  for  studying  single  or 
two-dcgrec-of-frccdom  flutter  problems. 
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